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1.  Introduction  and  summary. 

The  main  purpose  of  the  present  paper  is  to  give  a  rather  complete  tabulation  and 
discussion  of  properties  of  a  system  of  single  sampling  attribute  plans  dbtained 
by  minimising  average  costs  under  the  assumptions  that  costs  are  linear  in  p,  the 
fraction  defective,  and  that  the  distribution  of  lot  quality  is  a  double  binomial 
distribution. 

Starting  from  a  cost  function  containing  6  parameters  and  a  mixed  binomial  prior 

distribution  it  is  shown  how  the  average  costs  may  be  written  in  a  standard  form 

containing  only  two  parameters,  p  and  p  ,  besides  the  parameters  defining  the 

r  8 

prior  distribution.  The  one  parameter,  pr,  is  the  economic  break-even  quality  and 

depends  on  the  costs  of  acceptance  and  rejection  only,  whereas  the  second  parameter, 

p  ,  also  depends  on  the  costs  of  sampling  inspection  and  the  average  quality.  In  a 
8 

simple  and  practically  important  case  pr  and  p#  denote  the  costs  of  rejection  and 
the  cos ts  of  sampling  inspection,  respectively,  divided  by  the  costs  of  accepting 
a  defective  item. 

Specializing  the  prior  distribution  to  a  double  binomial  distribution  defined  by  the 
two  quality  levels  (p^,P2>  and  the  weights  (w^,w2),  Wj+  1,  it  will  be  seen  that 

the  optimum  sampling  plan  (n,c)  depends  on  the  6  parameters  (N,pr,pg,p1,p2,w2)  where 
N  denotes  lot  size.  It  cay  be  shown,  however,  that  the  weights  combine  with  the  p's 
in  such  a  way  that  only  5  (independent)  parameters  are  left. 

A  procedure  to  obtain  the  exact  solution  of  the  problem  has  been  developed  in  a 
previous  paper  and  this  has  been  used  for  computing  a  set  of  master  tables  in 
which  Pr«  pa«  0.01  tr.'l  0.10,  w2«  0.05,  ta^e  suitably  chosen  values  in 

relation  to  the  value  of  p^,  and  1  £  N  £  200,000. 
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la  the  remaining  part  of  Che  paper  Che  properties  of  the  optimum  plans  are  studied 
with  the  purpose  to  derive  simple  conversion  formulas  which  will  make  it  possible 
to  find  the  optimum  plan  for  an  arbitrary  set  of  parameters  from  a  plan  in  the 
master  table  with  a  "corresponding"  set  of  parameters.  The  main  tool  for  this 
investigation  is  the  asymptotic  expressions  for  the  acceptance  number  and  for  the 
sample  size,  viz. 

c  -  npQ+  a  +  o(l)  and  n  ■  —  (In  N  -  jlnln  N  +  lnX  +  j  ln<pQ)  +  o(l), 

<Po 

where  pQ  and  <p^  are  functions  of  (p^Pj)  only,  whereas  a  and  X.  depend  on  the  other 
parameters  also.  It  is  furthermore  shown  that  the  minimum  value  of  the  standardized 
costs  per  lot  asymtotically  equals  the  costs  of  sampling  inspection  plus  a  constant 
(depending  on  (p^,P2))  end  that  the  producer's  and  the  consumer's  risks  tend  to  zero 
inversely  proportional  to  lot  size.  Numerical  investigations  show  that  the  asymptotic 
expressions  give  good  approximations  to  the  optimum  plan  even  for  quite  small  values 
of  N. 


By  means  of  the  asymptotic  formulas  it  is  possible  to  find  out  how  (n,c)  vary  with 
the  individual  parameters.  One  of  the  most  Important  results  is  found  by  letting  all 
the  p's  tend  to  zero  which  leads  to  "the  proportionality  law":  The  optimum  sampling 


*  * 

plan  corresponding  to  (N, Xp  ,  Xp  ,  Xp.,Xp9,w9)  is  approximately  equal  to  (n  /\,c  )  where 
(n  ,c  )  is  the  plan  corresponding  to  (N  ,p  ,p  >P,,P*,w9)  with  N  »  NX. 

f  8  1  4  fc 

This  theorem  combined  with  other  similar  results  regarding  the  effect  of  varying  the 
individual  parameters  lead  to  two  general  conversion  formulas  stated  in  sections  8 
and  11.  A  summary  of  these  formulas  is  given  at  the  end  of  the  paper  in  connection 
with  the  tables. 


Efficiency  of  a  sampling  plan  is  defined  as  the  ratio  of  the  standardized  costs 
(loss)  of  the  optimum  plan  and  the  costs  of  the  plan  in  question.  Efficiency  is 
discussed  for  various  alternative  systems  and  the  efficiency  of  using  optimum  plans 
determined  fremwrong  values  of  the  parameters  is  studied. 

Finally  the  present  system  is  discussed  in  relation  to  other  systems  and  it  is 
pointed  out  that  from  an  economic  point  of  view  it  is  not  advisable  to  fix  the  con- 
or  the  producer's  risk.  If  one  wants  a  system  with  a  fixed  risk  then  the  risk 
should  be  fixed  to  50  per  cent  at  a  point  between  p^  and  Two  such  IQL  systems 
are  then  briefly  discussed. 


2 .  The  mode 1 . 

Several  authors  have  studies  economic  models,  mostly  linear,  for  the  determination 
of  single  sampling  inspection  plans  by  attributes,  see  for  Instance  [  1  ]  and  f  2  ] . 
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We  shall  here  stare  from  the  formulation  proposed  by  Guthrie  and  Johns  £  3  ]  and 
show  how  the  model  may  be  reduced  to  a  standard  form  as  previously  used  by  Hald 

O]. 

Let  N  and  n  denote  lot  size  and  sample  size  and  let  X  and  x  denote  number  of 
defectives  in  the  lot  and  the  sample,  respectively.  The  acceptance  number  is  denoted 
by  c. 


Let  the  costs  be 


and 


nSj+  xS£+  (N-n)Aj+  (X-x)A2  for  x  «  c  (1) 

nS^+  xSg+  (N-n)Rj+  (X-x)!^  for  x  >  c.  (2) 


The  interpretation  of  the  six  cost  parameters  depends  on  the  kind  of  Inspection 
envisaged,  i.e.  whether  inspection  is  a  consumers  receiving  inspection,  a  producers 
inspection  of  finished  goods,  or  "internal  inspection"  by  delivery  of  goods  from 
one  department  to  another  within  the  same  firm.  The  cost  parameters  may  have  quite 
different  values  when  considered  exclusively  froma  producer's  or  a  consumer's  point 
of  view  because  certain  costs  are  borne  primarily  by  one  of  the  parties  involved. 

The  values  of  the  cost  parameters  also  depend  on  whether  the  inspection  is  rectifying 
or  non-rectifying,  destructive  or  non-destructive.  In  the  following  the  two  cost 
expressions  are  discussed  and  a  few  examples  of  interpretation  are  given. 

Costs  associated  with  the  sample,  nS^d-  XS2,  for  brevity  called  "costs  of  sampling 
inspection"  consist  of  two  parts:  one  part,  nS^,  proportional  to  the  number  of 
items  in  the  sample  so  that  includes  sampling  and  testing  costs  per  item,  and 
another  part,  xS2>  proportional  to  the  number  of  defectives  in  the  sample,  i.e. 

S2  denotes  additional  costs  for  an  inspected  defective  item.  If  defective  items 
found  in  the  sample  are  repaired,  say,  then  S2  may  include  the  repair  costs  per 
item. 

"Costs  of  acceptance"  are  similarly  composed  of  a  part,  (N-n)A^,  proportional  to 
the  number  of  items  in  the  remainder  of  the  lot,  and  another  part,  (X-x)A2,  proportional 
to  the  number  of  defective  items  accepted.  Whereas  usually  will  be  zero  or 
negligible,  A2  will  often  be  considerable.  If  accepted  items  are  used  as  parts  in 
an  assembly  operation,  say,  A2  may  include  the  manufacturing  costs  (or  the  price) 
of  an  item,  the  costs  of  handling  the  defective  item  in  assembling  and  disassembling, 
and  the  damage  done  to  ocher  parts  used  in  the  assembly.  In  case  of  inspection  of 
finished  goods  A2  may  include  costs  of  repair,  service  and  guarantees  plus  loss  of 
good-will . 

"Costs  of  rejection"  consist  of  a  part,  (N-n)R^  proportional  to  the  number  of  items 
in  the  remainder  of  the  lot,  and  another  part,  (X-x)!^,  proportional  to  the  number  of 
defective  Items  rejected.  Rejection  is  here  taken  in  a  broad  sense  meaning  only  that 
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Che  lot  cannot  be  accepted  according  to  the  sampling  plan  used.  Rejection  may  there¬ 
fore  lead  to  sorting,  price  reduction,  scrapping,  or  salvaging.  If  rejection  means 
sorting,  say,  then  R^  Includes  sorting  costs  per  Item  and  R^  denotes  additional 
costs  for  defective  Items  found,  for  example  costs  of  repair  or  replacement. 

It  Is  obvious  that  from  a  practical  point  of  view  It  will  in  general  be  easiest 
to  obtain  information  on  the  values  of  the  cost  parameters  in  the  case  of  "internal 
Inspection". 

Denoting  the  hypergeometrlc  probability  by 


P(*|X)  -  $<£>/$ 


the  average  costs  for  lots  of  size  N  with  X  defectives  become 

n  c  n 

K(N,n,c,X)  -  E  (nS.+  xS_)p{x|X}  +  E  ((N-n)A.+(X-x)A0)  p(x|X}+  E  ((N-nJR.+fX-xjRpptxIXiO) 
x-0  1  2  x- 0  1  2  x-c+1  1 

Let  f^(X)  denote  the  (prior)  distribution  of  X,  i.e.  the  distribution  of  lot  quality. 

The  average  costs  then  become 

K(N,n,c)  -  E  K(N,n,c,X)fN(X).  (4) 

X 

As  shown  in  [  4  ]  this  expression  becomes  linear  in  N  for  the  important  class  of 
mixed  binomial  distributions,  i.e.  for 

fN(X)  -  /  <J)pV'X  dW(p)  (5) 

0 

where  W(p)  denotes  a  cumulative  distribution  function  (independent  of  N)  . 


From  (3)  -  (5)  we  find 


K(N,n,c)  -  /  K(N,n,c,p)dW(p) 

0 


where 


K(N,n,c,p)  -  n(Sj+  S2p)  +  (H-n)  ((Aj+  A2p)P(p)  +  (R^  R2p)Q(p))  , 


P(p)  -  B(c,n,p)  -  E  (")pXqn"X, 

x=0  * 


and  Q(p)  -  l-P(p)  . 

For  convenience  the  frequency  function  corresponding  to  W(p)  will  be  called  the 
distribution  of  the  process  average  or  the  distribution  of  p  as  distinct  from  fN(X) 


where 


* 
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which  gives  Che  distribution  of  X/N,  i.e.  the  distribution  of  lot  quality.  (The 
following  discussion  will  be  in  terms  of  p)  . 

Limiting  the  prior  distributions  to  mixed  binomials  (6)  shove  that  the  average 
costs  may  be  considered  as  an  average  of  the  cost  function  (7) .  which  is  a 
function  of  p,  with  respect  to  the  distribution  of  p.  It  should  be  noted  that 
this  result  is  valid  for  any  (N,n)  for  a  mixed  binomial  prior  distribution  and 
that  a  similar  result  holds  for  N  — >  «,  n  ->  »,  and  n/N  ~>  0,  for  any  prior 
distribution.  The  limit  theorems  derived  in  the  following  may  therefore  be 
applied  in  general. 

The  sampling  plans  discussed  are  obtained  by  minimizing  K(N,n,c)  according  to  (6) 
with  respect  to  (n,c)  for  given  cost  parameters  and  prior  distribution  and  they 
will  be  called  Bayesian  single  sampling  plans  or  optimum  plans. 

Starting  from  (7)  we  introduce  the  three  cost  functions 


k8(p>  - 

ka(P)  - 

and 

kr<P>  " 

defined  for  0  £  p  %  1.  We  shall  make  the 
functions : 


Sl+  S2P> 

(9) 

V  V’ 

(10) 

Rj+  R2P, 

(ID 

following  assumptions 

regarding  these 

1.  All  three  functions  are  non-negative  and  none  of  them  is  identical  zero. 

2.  k  (0)  <  k  (O)  and  k  (1)  >  k  (1),  from  which  follows  that  the  equation  k  (p)*k  (p) 

ci  it  dr  3  r 

has  the  solution 

pr=  (Rx-  r2),  0  <  pr  <  1,  (12) 


Pr  being  called  the  (economic)  break-even  quality. 

3.  kfi(p)  *  kra(p)  for  0  ^  p  5  1,  where 


km<P> 


J  ka<P> 

1  kr<P) 


for  p  5  pr 

for  p  >  p  . 

r 


(13) 


The  function  km(p)  gives  the  unavoidable  (minimum)  costs,  i.e.  the  costs  correspond¬ 
ing  to  the  situation  where  perfect  knowledge  of  quality  exists  without  costs  and 
all  lots  (processes)  are  classified  correctly  on  that  basis,  viz.  accepted  for 
p  5  pf  and  rejected  for  p  >  pf. 

Averages  over  the  prior  distribution  are  denoted  by  k  ,k  ,  etc.,  i.e. 

S  3. 


(14) 


-  6  - 


and 


ka-  /  ka(p)dW(p)  -  kfl(p)  «  Ax+  A2p, 

1  **r  1 

km“  /  km(p>dw(p)  "  /  k  (p)dW(p)  +  /  k  (p)dW(p). 
“  a  in  «  a  _  r 


(15) 


Costs  per  item  are  denoted  by  k,  costs  per  lot  by  the  corresponding  K,  i.e.  K  *  Nk. 

The  average  costs  for  the  three  cases  without  sampling  inspection,  i.e.  the  cases 
where 

(a)  all  lots  are  classified  correctly, 

(b)  all  lots  are  accepted,  and 

(c)  all  lots  are  rejected, 

then  become  k  ,  k  ,  and  k  ,  respectively.  These  cases  are  useful  "reference  cases" 
m  a  r 

since  sampling  inspection  is  justified  only  if  k  -  km<  min  [k^-  k^,  kf-  km), 
where  k  ■  K(N,n,c)/N. 

Case  (a)  will  usually  be  considered  as  the  basic  reference  case  and  average  costs 
for  other  cases  will  therefore  be  reduced  by  k  ,  since  k^  represents  the  average 
fixed  costs  per  item  which  will  be  incurred  irrespective  of  the  decision  made.  The 
cost  differences  , 


and 


ka“  km“  /  (ka(P)  *  kr(P))dW<P) 
Pr 


kr”  km=  /  <kr<P)  *  ka(p))dW(p) 


represent  average  decision  losses  in  case  (b)  and  (c)  respectively,  and  k  -  k 
represents  the  average  "loss"  by  Inspection. 

From  (6)  and  (15)  we  find 


and 


K  -  nk  +  (N-n)  /  (k  (p)P(p)  +  k  (p)Q(p))dM(p) 
s  0  a 


(16) 


leading  to 


K  -  nk  +  (N-n)  {  /  k  (p)dW(p)  +  /  k  (p)dW(p)} 
mm  i  a  r 

0  p 

rr 


*r  1 

K  -  K  -  n(k  -  k  )+(N-n)  (/  (k  (p) -k  (p))Q(p)dW(p)+  /  (k  (p)-k  (p))P(p)dW(p)  } 
m  sm  _ra  ar 

u  Pr 

^r  1 

*  n(k  -  k  )+(N-n)(A9-  R,)  {  /  (p  -  p)Q(p)dW(p)+  /  (p  -  P  ) P(p)dW(p)  }, 
sm  2  2  0  r  p  r 

the  two  terms  giving  the  average  costs  of  sampling  inspection  and  the  average 


(17) 
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Instead  of  minimizing  K  with  respect  to  (n,  c)  we  might  just  as  well  minimize  K-K^, 

(K-K  j/CiU-I!,,).  or  (K-K  )/(k  -k  ),  since  K  ,  A0-R0,  and  k  -k  are  independent  of  (n,c)» 

m  4  4  m  8m'  ID  4  4  8R1 

It  will  be  seen  from  (17)  that  it  is  practical  to  use  Av»R0  or  k  -k  as  "economic  unit".' 

t  t  s  m 

Defining  pr  1  pr 

P  “  /  pdW(p)  +  /  p  dW(p)  “  p  -  /  (p  -p)dW(p),  (18)  • 

m  0  p  0 


we  find  0  £ 

P  ~  P 

Km  *r 

and 

pr 

-  pm  * 
m 

(kr  -  km)/(A2  -  R,). 

(19) 

Defining  p 

8 

by  means 

of 

ps 

pm 

<ks  -  km>^A2  -  R2> 

(20) 

we  find 

ps 

-  Pr  « 

(ks  -  -  R2> 

(21) 

and 

P  * 
rs 

{(S1  -  Ax)+  (S2  -  R2)p)/(A2  -  R2). 

(22) 

Introducing 

R*(N,n 

,c) 

(K(N,n,c)  -  Km)/(A2  -  Rj) 

(23) 

we  find  the 

standard 

form 

pr 

1 

r  1 

r*  -  n(p  -  p  )  +  (N  -  n){  /  (p  -  p)Q(p)dW(p)  +  /  (p  -  p_)P(p)dW(p) )  (24) 

0  p 


containing  only  two  parameters  ?r  ans  pg  -  pm,  instead  of  the  six  cost  parameters  in 
the  original  model,  see  [4].  It  should  be  noted  that  p  -  p  depends  on  the  prior 

8  ID 

distribution  besides  on  the  cost  parameters. 

Consider  the  special  case  given  by  kfl(p)  ■  A^,  k^(p)  «  R^,  and  kg(p)  ■  which  is  a 

model  commonly  used  in  practice.  It  follows  that  p^  «  Rj/Aj  and  pg  *  ^1^2*  i.e.  pr  and 

p  are  the  costs  of  rejection  and  of  sampling  and  testing,  respectively,  measured  with 

the  cost  of  accepting  a  defective  item.  This  simple  interpretation  of  p  and  p  is  one 

r  8 

of  the  reasons  for  using  them  as  parameters. 

It  is  often  useful  to  discuss  the  problem  in  terms  of  the  simple  cost  functions 

k  (p)  =  p,  k  (p)  «  p  ,  and  k  (p)  *  p  ,  which  immediately  lead  to  the  form  (24).  The 
a  r  r  s  s 

corresponding  form  of  (7)  becomes 

KQ(p)  *  npg  +  (N  -  n)(pP(p)  +  PrQ(p)) 
from  which  the  general  form  may  be  found  as 

K(P)  a  <A2  -  R2)K0(p)  +  (nS2  +  (N  -  n)R2)p  +  (NA^  -  n(S2-  R2)p). 

A  sketch  of  the  cost  functions  for  a  typical  case  has  been  given  in  Fig.  1,  which 

is  bn.9od  on  thr  data  in  ncrtion  13. 


Fig.  1 

Example  of  cost  functions. 
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For  some  purposes  it  is  useful  to  use  kg-  km  as  economic  unit  instead  of  A^-  R2 . 
Putting 


R(N,n,c)  -  (K(N,n,c)  km). 


i  .e. 


R  -  R*/(ps- 


> 


we  find 

N_  Pr  * 

R  s  n  +  ~~~p  (  /  (Pr-  p)Q(p)dW(p)  +  /  (p  -  Pr)P(p)dW(p)  },  (25) 

*s  0  p 

rr 

the  two  terms  again  giving  the  costs  of  sampling  inspection  and  the  average  decision 
losses,  respectively,  but  here  using  the  average  costs  of  sampling  inspection  (minus 
km)  per  item  in  the  sample  as  economic  unit. 


In  the  next  section  we  shall  discuss  the  determination  of  (n,c)  for  a  double 
binomial  distribution  as  prior  distribution.  This  means  that  p  is  a  random  variable 
taking  on  only  two  values,  p^<  p^<  p 2,  with  probabilities  w^  and  w^  1  -  w^, 
respectively.  From  (25)  we  then  find 

R  -  n  +  (N  -  n)(71Q(p1)  +  72P(p2))  (26) 


where 


V  'V  pr'wi/(ps-  Pm)  =  lka(pi>’  kr<pi>  I  V°V  V'  l'1'2'  (27) 


piwi+ 


P  wo 
r  2 


(28) 


i.e.  R  depends  on  four  parameters  only,  viz.  P^,P2,7j^72. 

The  correspondingly  standardized  costs  for  the  cases  of  acceptance  and  rejection 
without  inspection  are 

V  N<v  V'V  v  ■  **2  <29> 


and 

Rr"  N<V  km>  *  "’’l*  <*» 

I 

These  results  may  also  be  obtained  from  (26)  for  n  =  0  by  setting  P(p)  =  1  and  0, 
respectively. 

If  acceptance  without  inspection  is  cheaper  than  rejection  without  inspection,  i.e. 

k  <  k  we  find  p  <  p  and  y0<  7.  . 
nr  r  z  1 


/ 
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In  the  special  case  k  ■  k  we  have  p  ■  p  and  7  *  1  so  that  the  model  contains  only 

St  S  IT  1 

three  parameters. 


It  should  be  noted  that 


and 


1 


1 


(31) 


(32) 


3.  The  exact  solution  and  the  tables. 

In  a  previous  paper  [  4  ]  we  have  proved  the  following  theorem: 

For  a  double  binomial  (prior)  distribution  of  lot  quality  Riven  by  the  parameters 
(Pl,p2,w2)  and  for  linear  cost  functions  (1)  and  (2)  the  Bayesian  single  sampling 
plan  may  be  found  by  minimizing  R  (N,n,c),  see  (26),  with  respect  to  (n,c)  .  The 
solution  satisfies  the  two  inequalities 


and 

where 


and 


a+|3cgn<or  +  p(c  +  1) 


F(n  -  l,c)  S  N  <  F(n,c) 


log 


W2(P2~  Pr)  / 
wl(pr"  pl*  ' 


log  —  =  log 


72  /  ql 
7i  / 108  q2  ' 


3  -  log 


p2ql  /  .  ql 
q2pl  '  q2  ’ 


F(n,c)  =  n 


P  -  P  +  E  w  (p-  p  )B(c,n,p  ) 
s  r  ^  1  t  1  1 

+  1  +  £  w  (p  -  Pr>Pib(c^n#Pi) 


(33) 

(34) 

(35) 

(36) 

(37) 


tr 


For  two  plans  (n^Cj)  and  (n2,c2)  ,  c1<  c2  say,  satisfying  C33)  and  having  overlapping 
N-intervals  according  to  (34)  R(N,n1,c1)^  R(N,n2,c2)  for  N  ^  N12  where 


N 


12 


(Ps-  PrHn2*  ni>  +  n27(n2'c2^*nl7^nl,cl^ 
7(n2,c2)  -  7(1^,^) 


(38) 


and 

7(n,c)  =  £  wi(pr-pi)B(c,n,pi).  (39) 

i 

In  [  4  J  the  theorem  was  derived  as  a  special  case  of  a  more  general  one.  We  shall 
here  derive  the  theorem  directly  from  (26)  using  the  same  method  as  in  [  4  ] . 

Values  of  (n,c)  minimizing  R  must  satisfy  the  two  inequalities 
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£cR(N,c-l,n)  *  0  <  £cR(N,c,n),  O  *  c  *  n,  (40) 

and 

AnR(N,c,n-l)  S  0  <  £>nR(N,c,n),  c  *  n  ;?  N,  (41) 

A  denoting  the  usual  forward  difference  operator. 

Noting  that  AcB(c,n,p)  -  b(c  +  l,n,p)  and  AQB(c,n,p)  *  -pb(c,n,p)  we  find  from  (26) 

A£R(N,n,o)  -  (N-n)  (-7^(0+!, n,pL)  +  72b(c+l,n,p2) }  (42) 

and 

AnR(N,n,c)  =  l-  {71Q(P1)+72p(P2))  +  (N-n-1 ) (71P1b (c, n, px )-72P2b (c,n, p2 )  ). (43) 

Inserting  these  expressions  into  (40)  and  (41)  and  solving  for  n  and  N,  respectively, 
immediately  leads  to  (33)  and  (34).  From  R(N,n^,c^)  =  R(N,  n2,c2)  we  next  determine 
^12  ky  solving  for  N. 

A  sketch  of  R  as  function  of  n  and  c  for  fixed  N  has  been  given  in  Fig.  2  for  a 
typical  case. 

The  economic  interpretation  of  (40)  and  (42)  is  the  following:  For  given  n  the  optimum 
value  of  c  is  determined  such  that  a  change  of  c,  an  increase  by  1  say,  will  give 
nearly  no  change  of  the  total  decision  loss,  since  the  loss  due  to  the  increased 
consumers  risk  is  nearly  balanced  by  the  gain  due  to  the  smaller  producer's  risk. 

Similarly  the  interpretation  of  (41)  and  (43)  is  that  for  given  c  the  optimum  value 
of  n  is  determined  such  that  a  change  of  n,  an  increase  by  1  say,  will  give  nearly  no 
change  of  total  costs,  since  the  increase  of  sampling  inspection  costs  by  1  minus 
the  average  decision  loss  for  one  item  is  nearly  balanced  by  the  decrease  in  decision 
losses  for  the  remainder  of.  the  lot. 

Tabulation  of  optimum  plans  may  be  carried  out  by  starting  from  the  smallest  value  of 
c  giving  a  positive  n(n  >  c)  according  to  (33),  i.e.  c^  =  [-  a/([3-l)J,  [  ]  denoting 

"the  integer  part  of".  For  consecutive  values  of  c,  n-  and  N-intcrvals  are  computed 
from  (33)  and  (34)  and  in  case  of  overlapping  N-intervals  costs  are  compared  by 
means  of  (38).  A  detailed  example  may  be  found  in  [4],  The  tables  have  been  computed 
by  thi3  method  on  an  electronic  computer. 

The  sampling  plans  have  been  tabulated  for  two  "quality  levels",  viz.  p  =  p  =  0„0J 

r  s 

and  0.10,  for  one  value  of  the  weight  function  w2  =  0.05,  for  8  values  of  p^/p^,  and 
for  10,  respectively  5,  values  of  p2/pr,  giving  a  total  of  120  tables.  Each  table 
gives  (n,c)  as  function  of  N  for  N  <  200,000. 


-  Zi- 
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For  =  0.01  the  search  for  optimum  plans  has  been  limited  to  values  of  n  which  ~,-c 
multiples  of  5, 

These  tables  will  be  referred  to  as  "master  tables"  since  optimum  plans  for  other 
values  of  the  parameters  may  easily  be  found  from  the  tabulated  ones  by  means  of 
conversion  formulas  developed  in  the  following  sections. 

The  exact  solution  has  been  modified  in  one  respect.  For  a  given  value  of  c 

the  first  and  last  N-interval  may  be  rather  short  as  compared  to  the  other  intervals. 

As  an  example  consider  the  following  section  of  the  original  table  for  p  «*  p  ~  0,010, 
p^*  0.006,  p2  =  0,020,  w2  =  0.05: 


II  n  c  £N 


4010  - 

4370 

165 

3 

360 

4370  - 

4420 

170 

3 

50 

4420  - 

4430 

240 

4 

10 

4430  - 

4920 

245 

4 

490 

4920  - 

5570 

250 

4 

650 

5570  - 

5590 

255 

4 

20 

5590  - 

5610 

325 

5 

20 

5610  - 

6250 

330 

5 

640 

The  example  shown  is  an  extreme  one  with  small  intervals  occurring  at  the  beginning  as 
Well  as  at  the  end  of  each  section  of  the  table.  It  is  naturally  without  any  interest 
to  use  the  sampling  plan  (240,4)  for  4420  <  N  <  4430  and  then  change  to  (245,4)  for 
4430  <  N  <  4920.  To  eliminate  such  small  intervals  from  the  final  table  it  was  decided 
to  discard  the  first  and  the  last  sampling  plan  for  a  given  c  if  the  length  of  the 
Corresponding  N-intcrval  was  less  than  1/5  of  the  length  of  the  neighbouring  interval. 
In  such  cases  the  value  of  N  according  to  (38)  was  computed  for  the  new  neighbouring 
plans,  (165,3)  and  (245,4)  say,  to  find  the  optimum  N-intervals  for  the  remaining 
plans.  The  result  of  this  procedure  is  in  most  cases  practically  equal  to  incorporating 
the  small  N-intervals  into  the  larger  neighbouring  intervals,  for  example  using  (245,4) 
for  4420  <  N  <  4920. 

To  save  space  every  second  N-interval  for  a  given  value  of  c  has  been  omitted  because 
mrresponding  sampling  plans  may  be  found  by  addi  ng  1  (p  =  0.10)  and  5  (p  =  0.01), 
lyspectivel y .  to  n  for  the  preceding  interval. 

▼Slues  of  N  have  been  rounded  to  3  significant  figures  and  tabulation  has  been  stopped 
«t  N  =»  200,000. 

AM  mentioned  above  the  tables  were  designed  as  master  tables  from  which  optimum  plans 
nay  be  derived  for  other  values  of  the  parameters  and  for  this  reason  it  was  decided 
t©  tabulate  the  complete  solution  with  respect  to  N  to  make  interpolation  superfluous. 


t 


The  user  of  the  tables  in  practice  may  easily  derive  a  simplified  set  of  tables 
from  the  given  ones,  either  by  using  a  set  of  fixed  N-intervals,  or  a  set  of 
fixed  N- arguments.  An  example  has  been  given  in  the  following  table. 

Single  Sampling  Tables  for  100pr  r  100pg  -  1.0,  100p^  ~  0.5,  and  W2  -  0.05. 


100p2 

1. 

5 

2. 

0 

3. 

0 

4.0 

5.0 

6.0 

-  - ' 

7.0 

N 

« 

n 

c 

n 

c 

n 

c 

n 

1 

c 

n 

c 

n 

i 

c  j 

n 

c 

20 

j 

Accept 

r 

0 

30 

Accept 

5 

0  1 

1 

5 

0 

50 

5 

0 

10 

0  i 

10 

0 

70 

5 

0 

10 

0 

15 

0 

100 

Accept 

10 

0 

15 

0 

15 

0 

200 

10 

0 

20 

0 

20 

0 

25 

0 

300 

15 

0 

50 

1 

50 

1 

50 

1 

500 

Accept 

55 

1 

60 

1 

60 

1 

55 

1 

700 

45 

1 

65 

1 

65 

1 

65 

1 

60 

1 

1CXX> 

55 

1 

110 

2 

105 

2 

100 

2 

90 

2 

2000 

125 

2 

170 

3 

155 

3 

140 

3 

lOO 

2 

3000 

Aeoept 

195 

3 

180 

3 

160 

3 

145 

3 

135 

3 

5000 

185 

3 

265 

4 

235 

4 

210 

4 

185 

4 

140 

3 

7000 

275 

4 

330 

5 

290 

5 

215 

4 

190 

4 

170 

4 

10000 

450 

6 

400 

6 

295 

5 

260 

5 

195 

4 

175 

4 

20000 

I  Accept 

640 

8 

475 

7 

355 

6 

310 

6 

240 

5 

215 

5 

30000 

755 

9 

735 

9 

545 

8 

410 

7 

315 

6 

280 

6 

220 

5 

50000 

1080 

12 

920 

11 

620 

9 

470 

8 

365 

7 

285 

6 

255 

6 

70000 

1300 

14 

1095 

13 

690 

10 

475 

8 

410 

8 

325 

7 

260 

6 

100000 

1520 

16 

1190 

1 

14 

760 

11 

530 

9 

415 

8 

330 

7 

290 

7 

The 'haturaf' parameters  of  the  model  are  (Pi/P2*w2^  which  characterize  the  prior 
distribution,  and  (p  ,p  ),  which  depend  on  the  costs.  The  tables  and  the  properties 
of  the  solution  will  be  discussed  in  terms  of  these  parameters  on  basis  of  the 


results  in  the  next  section.  However,  one  property  may  be  stated  immediately  from 

the  observation  that  the  solution  depends  on  four  parameters  only,  viz.  (p^,P2>7^>72^* 

The  three  parameters  (p  ,p  ,w«)  may  therefore  in  respect  to  the  solution  be 

considered  as  functionally  related,  i.e.  combinations  of  (p  ,p  ,w_)  Riving  the  same 

1  r  s  / 

viU  lead  to  the  same  sampling  plan. 

From 


we  find 


From 


1 


V.  -  pm> 


<p2-pr)w2 
<pr  pl)vl 


(P< 


P,)/(l  + 


?2W1 

71W2 


). 


(Pr  - 


(44) 


and 

we  f ind 


m 


Ps  '  P1 


P1  '  w2<pr 


pl> 


w 


1 


<pr  -  Pl><—  +«2>. 


(45) 


*2*  *0  ^  XS.  -IN  Ol 

^  ^  *-  •  x  rK  rv 


K)  Ch  Ot 


^  K  % 
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These  formulas  show  how  and  pg  depend  on  ^or  8iven  (Pi# use 


them  in  connection  with  the  master  tables  we  put  p 


p  and  w. 
ro  2 


0.05a.  which  leads 


to 


207 


pr<x>  -  pl0+  <p2o-  p1o)  /a  -  r2o*  —  > 


where 


P2o'  pro  p2~  1  , 

2°  1»<P„-  Plo)  '  19(1  -  0l) 


the  index  o  denoting  an  argument  in  the  master  table,  Pro“  0.01  or  0.10,  p£0/Pro 
Dividing  by  p  gives 


pr(  )/pro“  cl+  (o2'  pl)  '  (1  *  ?,2o+ 
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2o 


)  * 


(46) 


which  has  been  tabjlated  in  the  appendix. 

The  field  of  application  of  the  master  tables  may  therefore  be  considerably  enlarged 
by  making  use  of  the  following  rule: 

The  optimum  sampling  plan  for  (M,  pfQ, p^,  P2o,w2*  Pro*  ^so*  l-3  the  same  as  the 

Ela.ri.ior  <N>Prof <w2'pl'p2)  ,plo,p2o,W2)  * 

Consider  for  example  the  case  with  p  «  p  -  0.010.  p.-  0.006,  p  *  0.040,  and 

ro  so  i  Z 

W2"  0.05  for  which  the  optimum  plans  have  been  given  in  the  master  table.  The  same 
plans  are  also  optimum  for  w  *  0.20,  say,  and  p  *  p  =  0.019,  p*  0.006,  and  p  =  0.040 
which  may  be  seen  by  interpolation  in  the  table  of  f(«2,p1,P2)  for  p^  0.6  and 
4.0. 


4.  The  asymptotic  solution. 

In  this  section  we  shall  give  a  somewhat  simpler  and  more  direct  proof  of  the 
asymptotic  results  found  by  Guthrie  and  Johns  [3]  and  by  Ha  Id  {”  4  ] ,  and  further¬ 
more  carry  the  asymptotic  expansion  so  far  that  we  get  a  useful  approximation  to 
the  exact  solution  also  for  small  values  of  c. 

The  proof  is  based  on  the  following  lemma  which  is  a  special  case  of  a  theorem 
proved  by  Blackwell  and  Hodges  [5]: 

For  c/n  =  h  =  p  +  t,  p  being  a  constant  and  e  — >  0  for  n  — >  00,  we  have 
o  o 

P(P)  ■  h  -daL  c-nx>(h,p)(1  +  0(l^))  £or  p  <  P) 

]/  2  imp  q  p  -  p 

’  00  o 


(47) 


where 


<p<h,p)  -  h  In  (1  -  h) In 


(48) 


For  pQ>  p  the  tome  expression  is  valid  for  Q(p)  if  only  p  -  pQ  is  replaced  by  pQ-  p. 


Writing 


where 


<P<h>P)  -  <p(po,p)  +  cy’ivo,p)  +  0(e  ) 


<p'(P0>P) 


po<* 

V 


(49) 


(50) 


we  find  from  (26)  and  (47)  the  asymptotic  expression 

<U 


R  -  n  +  (N-n) 


/  2rmp 


2 

I 

o«o  1-1 


hpt 


Po"  pi 


-ncp(p  ,  p. )  -nc c p’  (p  ,  p. ) 
e  °  1  °  1  (1  +  0(#)) 


(51) 


on  the  assumption  that  p.<  p  <  p_.  (As  will  be  shown  later  e  -  0(l/n),  and  we  may 

ty  A  O 

therefore  disregard  ne  )  .  We  shall  first  determine  the  value  of  h  *  pQ+  e  which 
minimize  R  for  given  n  and  next  determine  the  value  of  n  giving  the  absolute  minimum 
by  treating  R  as  a  differentiable  function  of  h  and  n. 


The  essential  feature  of  (51)  is  that  the  two  binomial  risks,  Q(p^)  and  P(p2),  have 
been  expressed  as  functions  tending  exponentially  to  zero  for  n  — >  «. 

As  explained  in  [4  ]  the  optimum  plan  must  have  the  property  that  R/N  — >  O  for 
N  ->  »,  n  “>  oo,  and  n/N  — >  0.  It  follows  that  pQ  must  satisfy  the  Inequality 
p^  <  pQ<  p^  because  otherwise  R/N  would  not  tend  to  zero  but  to  y^  or  y^. 

We  shall  state  the  theorem  to  be  proved  for  the  double  binomial  distribution  only, 

but  it  is  valid  for  a  more  general  class  of  distributions,  viz.  for  a  distribution 

having  probability  density  w(p)  »  0  for  p^  <  p  <  p2,  w(p^)  *  Wj>  0,  w(p^)  •  w2>  0, 

Wj+  w„  £  1,  and 
A  ^  * 

pl  1 

/  dW(p)  +  /  dW(p)  -  l-w.-w, 

0  * 
p2 

#  * 

for  0  £  Pj<  j ;  and  p2<  r  .,£  1,  which  means  that  the  probability  distribution  may  be 
arbitrary  outside  the  interval  p^<  p  <  p2«  The  result  of  such  a  generalization  will 
only  be  to  add  a  term  to  (51)  of  form 


N-n 


r 


ps  Pra  ^2»mpoqo 


(P  ‘  P)P 
/  — S -  e 

I  po’  p 


“n<P(h>P)dW(p), 


(I  denoting  the  intervals  (0  £  p  £  p^)  and  (p2  £  p  £  1))  which  obviously  is  0(e’n) 


times  the  last  term  of  (51)  since  4>(h,p)  >  <p(h,p^)  for  p  <  P^  and  cp(h,p)  >  'p(h> 


for  p  >  Pj. 


Because  of  the  factor  pf-  p  in  the  cost  function  we  might  also  have  assumed  that 
w(pr>  >  0  without  altering  the  result. 

It  is  reasonable  to  assume  that  the  two  exponential  terms  in  (51)  tend  to  zero  with 


the  same  speed,  l.e. 

that  p  is  determined 
ro 

from 

f(P0^Pl) 

"  9(P0>P2) 

which  gives 

°  q2  / 

;  p2qi  i 

q2P!  0 

(52) 

and 

i  P° 

V  poln  57  ■ 

q 

4-  q  ln  ~  ,  i  *  1  or  2 . 

Ho  q.  * 

(53) 

Under  this  assumption  we  shall  determine  e  by  minimization  of  (51)  .  The  part  of  R 
depending  on  e  is 

...  .  Vi  -«*’<p0>pt>  . 

f(e)  *«  I  -  e 


From  f *  (e)  *  O  we  find 


1  iv  pi 


2  7  p  .I 

1  -  e  =0 


i-1  p  -  p. 

1  o  1 


(54) 


where  -  according  to  (50)  - 


<P[  “  ln 


Poqi 

Vi 


(55) 


Solving  for  a  *  ne  we  find 


where 


a5* 

o 


ln 


ylpl(p2~  Po)vl 
y2P2(po"  Pl><-f02> 


do  =  91  ‘  ^2  =  ln 


p2ql 

q2Pl 


(56) 


(57) 


We  thus  have  the  result  that  c  =  npQ+  a  +o(l)  in  accordance  with  what  could  be  expected 
from  (33) . 


Inserting  these  results  into  (51)  we  find 


-n«p 

R  *  n  +  (N-n)  e 


with 


2  7iPl  -ap' 

I  -  e 


V  2rtp  q  i*=l  Ip  -  p. 
*oMo  1  O  1 


(58) 


(59) 
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To  prove  (Indirectly)  that  h  ■  pQ+  e  minimizes  R  let  us  assume  that  h  ■  pQ+  e, 
given  by  (52)  and  (56),  does  not  minimize  R  but  that  min  R  is  obtained  for 
h  »  ho+  cQt  hQ4  Pc  and  co™>  ®®not^n8  the  part  of  R  depending  on  h  by  g(h)  we 

find  for  sufficiently  large  n  and  for  hQ<  p  ,  say ,  that 


-n<p(h  ,p  ) 

g(h0)  -  Xj(ho)  e  0  1  (1  +  0(e  ")) 


since  cp(h  ,p.)  >  <p(h  ,p.)  for  h  <  p  .  However.  g(h  )  cannot  be  min  g(h)  since  cp(h  ,p.) 

Ol  O  1  o  o  o  O  JL 

<  <p(p  »p,)>  i.e.  we  have  reached  a  contradiction  by  assuming  h  4  P  • 

OX  o 1  o 


From  dR/dn  *  0  we  find 


1  -  (N-n)  ft  e  %o+  y  -  ft  e 


(60) 


ln(N-n)  «  cpQ n  +  j  In  n  -  ln(A.q>o)  +  o(l) 


(61) 


From  (58)  and  (60)  we  also  have  that 


min  R  =  n  +  -  +  o(l) 

(n,c)  ^o 


(62) 


where  n  may  be  determined  by  inversion  of  (61),  i.e. 

n  ■  ~  (In  N  -  ^  In  In  N  +  In  \  +  r  In  if  )  +  o(l)  . 

(p  t  to 

o 

We  have  thus  found  that  asymptotically  c  is  a  linear  function  of  n  and  n  is  pro- 
portional  to  In  N-  j  In  In  N  plus  a  constant.  Furthermore  it  follows  from  (62)  that 
the  average  decision  loss  per  lot  tends  to  a  constant  l/tpo  so  that  for  large  lots 
decision  losses  divided  by  sampling  inspection  costs  tend  to  zero. 

To  investigate  the  two  risks  asymptotically  we  find  from  (54) 


71P1<P1  "a^l  72p2^”cpP  ’atf2 


V  pl 


p2-  Po 


so  that  (59)  gives 


%  7iPi8o  -a“’i 

x  <p°-  pi)(-’P  e 

’  ro  o 


which  together  with  (60)  may  be  used  to  reduce 


ru  ^  1  qoPl 

Q(Pl  ’  p°"  pi 

’  o^o 


■acpj  -n<pc 


1 

N-n  * 


(63) 


Similarly  we  have 


<,(Pl)  "  9a/,a‘ 

O  1  O 


Kp,) 


V26i  N*n 


so  that 


p(p2)/Q(p1)  *  71cPi/r2(“‘Pp. 


(64) 


(65) 


We  have  thus  proved  the  following  theorem; 

Asymptotically  the  optimum  sampling  plan  is  given  by 


and 


which  lead  to 


n 


min  R 


Q(PX) 

and 

P(P2) 


c  -  npQ  +  a  +  o(l) 

-  i-(ln  N-|  In  In  N  +  In  X+|  In  tp  )  +  o(l) 

<p  i  i  o 

*  ~  (In  N  In  In  N  +  In  X  +  j  In  +  1)4-  o(l), 

CPq 


m*2 


cp  7.8' 
o  1  o 


N-n  +  0(N> 


<P 1 


cp  7o$ 

o  2  o 


1  ^  A\ 

N-n  +  °(N*' 


(66) 


(67) 

(68) 


It  will  be  noted  that  pq  and  depend  on  (P^>P2)  only,  i.e.  they  are  independent 
of  the  cost  parameters  and  of  w^. 

The  asymptotic  solution  supplements  the  exact  one  in  several  respects.  Since  the 
optimum  plan  is  a  function  of  5  parameters  (N, p^,Pj >7^,72 )  a  complete  tabulation 
is  rather  hopeless  even  if  a  program  has  been  worked  out  for  an  electronic  computer. 
Furthermore  the  properties  of  the  exact  solution  are  not  easily  to  be  found  from 
the  procedure  by  which  the  solution  is  obtained.  The  advantages  of  the  asymptotic 
solution  are  that 


(1)  it  clearly  shows  how  the  optimum  plan  and  various  derived  quantities  depend 
on  the  parameters t 

(2)  it  may  be  used  as  starting  point  for  developing  approximations  which  are  valid 
also  for  small  N, 

(3)  it  may  be  used  for  developing  interpolation  and  extrapolation  formulas  in 
connection  with  "master  tables”  of  the  exact  solution,  and 

(4)  it  shows  the  sensitivity  of  the  solution  with  respect  to  changes  of  the 


parameters. 
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These  aspects  of  the  solution  will  be  discussed  in  the  following  sections. 

4 


5.  Comparison  of  exact  and  approximate  solution. 

Looking  at  the  relation  between  n  and  c  in  the  tables  it  will  be  seen  that  the 
optimum  values  of  n  for  a  given  value  of  c  tend  to  cluster  around 


nc  =  a  +  (3  (c  +  j) 


(69) 


as  might  be  expected  from  (33).  Comparing  with  the  asymptotic  result  c  =  npQ  +  a, 
Po  =  1/(3  and  a  being  defined  by  (56),  agreement  between  the  two  expressions  would 
require  that 


p2(Po-Pi>(V2>\  //  p2qx\ 

Pl^-PoS 


1 

2 


It  can  be  proved  that  the  ratio  on  the  left  hand  side  above  is  positive  and  less 
that  1.  Numerical  investigations  show  that  in  typical  coses  in  practice  the  ratio  does 
not  deviate  much  from  1/2.  As  examples  consider  the  following  results; 


100?! 

100p2 

p2/pl 

Ratio 

0.2 

4.0 

20 

0.  528 

0.2 

2.0 

10 

0.517 

0.6 

4.0 

6.7 

0.512 

0.  6 

2.0 

3.3 

0.505 

The  ratio  depends  primarily  on  p2/p^  cn(^  practically  the  same  results  will  be  found 
for  values  offp^p^  which  are  10  times  as  large  or  l/lo  of  the  values  considered. 

We  shall  therefore  in  the  following  use  the  simpler  expression  (69)  instead  of 
c  “  npQ  +  a  as  the  starting  point  for  finding  n  from  c  or  reversely. 

The  asymptotic  formulas  may  be  used  in  two  ways 

(1)  Starting  from  c  we  may  determine  the  corresponding  N-interval  and  within 

l  that  the  relation  between  n  and  N. 

(2)  Starting  from  N  we  may  determine  the  corresponding  n  and  from  n  determine  c. 

The  first  method  is  useful  for  making  a  systematic  tabulation  of  sampling  plans 
whereas  the  second  is  suitable  for  computing  "isolated"  plans  for  a  given  N. 

Starting  from  an  integer  value  of  c  we  flrut  find  nc  from  (69)  and  the  corresponding 

N  from  (61).  Similarly  we  find  N  _  .  and  N  _  c,  being  the  lower  and  upper  limit 
c  J  c-0. 5  c+0.  5 

for  N  having  c  as  optimum  acceptance  number. 
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In  the  asymptotic  solution  we  have  disregarded  the  discreteness  of  c  and  n.  We  may, 
however,  afterwards  try  to  take  the  effect  of  the  discreteness  of  c  into  account 
by  investigating  the  relationship  between  n  and  N  for  given  (integer)  value  of  c. 
From  dR(N,n,c)/dn  =  0  it  can  be  found  that  n  is  approximately  a  linear  function  of 
log  N  with  slope  -1/log  *)  .  Within  the  interval  (Nc_q  5>Nc+q  5)  wc  133 y  therefore 
determine  n  from  the  approximate  formula 


n  =  nc-  (log  N  -  log  Nc)/log  q2,  NC_0  5<N  <  Nc+0.5' 


(70) 


which  for  small  p2  and  small  intervals  nay  be  replaced  by 


n  -  n  +  (N  -  N  )/N  p0,  N  _  ,<  N  <  N 

c  '  c  cr2'  c-0.5  c+0.5 


(71) 


It  follows  Lhat  the  values  of  n  belong  to  the  interval 

■V  P('o+  2T")/2l>2- 
c 

For  applications  in  practice  we  give  the  formula  corresponding  to  (61)  with 
logarithms  to  base  10.  i.e. 


where 


log  (Nc-  nc)  =  cf)nc+  -  log  nc  + 


(72) 


c?  =  p  log  -  +  q  log  , 

o  p.  ^ o  q. 

=  -log  (Xtpo) , 


1  or  2, 


(73) 


(74) 


<p(a+f>  '  q  2 

Xf  =  10  2  ■■■*  ~  Jr-2-  S 

o  log  e  \l 2  tp  . 

0  ^  o  1”1 


7ipi 


p  -  p. 
o  1 


(75) 


-a  having  been  replaced  by  ^  j  ill  Xq>o . 

In  the  following  wc  shall  make  much  use  of  (72)  with  N^-  n£  replaced  by  N  which 
only  means  that  wc  disregard  terms  of  order  nc/Nc  and  less. 

The  approximation  obtained  by  using  (69),  (70),  and  (72)  is  usually  very  good  even 
for  quite  small  values  of  c.  Normally  the  approximate  value  of  c  will  deviate  at 
most  1  from  the  correct  value.  The  approximation  depends  essentially  on  p^/p  , 
being  good  for  large  values  of  p2/p^  and  poorer  for  snail  values.  Two  examples 
for  p2/pj  =  6.7  ancl  3.3,  respectively,  will  show  the  results  obtained  for  a  typical 
good  and  poor  case.  Table  1  and  Fig.  3  show  that  the  approximate  and  the  exact 
solution  are  practically  identical  in  the  first 


This  result  is  due  to  Mrs.  K.  West  Andersen. 
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case  whereas  the  approximate  solution  in  the  second  case  often  will  lead  to  a  value 
of  c  being  1  too  large  and  a  corresponding  value  of  n. 

Table  1. 

Comparisons  of  exact  and  approximate  sampling  plans  computed  from  (69)  ,(70),  and  (72). 


Pr=  pn=  0.010,  p^=  0.006,  p2=  0.040,  w^  0.05. 
a  =  -26.7,  p  =  55.509,  *  0.0034156,  5  =  1.5499,  -1/log  q2  «=  56.405. 


Approximation 

Exact 

c 

n 

c 

n 

Nc±0 . 5 

n 

N 

1 

57 

43-66 

269-714 

45-65 

280-714 

2 

112 

104-120 

715-1400 

105-120 

715-1420 

4 

223 

216-230 

2490-4300 

220-230 

2550-4390 

6 

334 

320-340 

7190-12000 

330-340 

7390-12300 

8 

445 

439-451 

19700-32300 

440-450 

20200-33000 

10 

556 

550-562 

52400-85300 

550-560 

53600-87000 

12 

667 

661-673 

137000-200000 

665-670 

140000-200000 

pr=  Pg=  0.010 

,  px=  0.006,  p2= 

:  0.020,  w2- 

0.05. 

a  «  -143.0,  p  =  85.879,  <p  =  0.0009088,  8  =  2.0785,  -1/log  q^  113.97. 


Approximation 

Exact 

c 

n 

c 

n 

Ws 

n 

N 

2 

72 

44-83 

715-1750 

- 

- 

4 

243 

234-251 

2790-3970 

245-250 

4420-5590 

6 

415 

408-422 

5410-7150 

415-420 

7100-8980 

8 

587 

581-593 

9270-11900 

585-595 

11300-14200 

10 

759 

753-765 

15100-19000 

755-765 

17700-22000 

12 

Ml 

925-936 

23800-29600 

930-935 

27300-33700 

14 

1102 

1097-1107 

36800-45700 

1100-1105 

41500-51100 

16 

|  1274 

1269-1279 

56500-69700 

1270-1280 

62800-77200 

18 

1446 

1441-1451 

86000-106000 

1445-1450 

94600-116000 

20 

1618 

1613-1623 

130000-159000 

1615-1620 

142000-173000 

Comparisons  of  exact  and  approximate  sampling  plans. 


It  io  essential  for  the  efficiency  of  the  approximation  to  use  Che  right  relation 
between  n  and  c,  see  tie  discussion  in  section  12,  and  it  is  therefore  fortunate 
that  this  relation  is  a  simple  linear  one. 


The  approximation  formula  breaks  down  for  values  of  N  for  which  the  cheapest  solution 
is  acceptance  without  inspection  (or  rejection  without  inspection).  As  will  be  seen 
from  Table  1  the  approximation  formula  may  in  such  cases  lead  to  a  sampling  plan  even 
if  no  optimum  plan  exists.  The  difference  in  costs  by  using  such  a  plan  instead  of 
accepting  without  sampling  inspection  will,  however,  normally  be  small. 


Turning  to  the  inverse  formula  (67)  numerical  investigations  show  that  the  results 
are  not  as  accurate  as  those  found  from  (61)  .  Taking  one  more  term  in  the  inversion  of 
(61)  and  changing  to  logarithms  with  base  10  we  find 


where 


V 


;  uog  n 


<2  log  log  N  +  d)  (1  ------ 


)] 


(76) 


d  =  -  log  \cpQ  -  j  log  cp 


6  -  —  log  cp. 


(77) 


The  exact  Inversion  leads  to  the  correction  term  (log  e)^2  log  N)=  0.22/log  N  which, 
however,  on  the  basis  of  numerical  investigations  has  been  replaced  by  1/p  log  N)  . 

If  (76)  is  to  be  used  extensively  it  pays  to  tabulate 

g(N)  -  log  N  -  j(l  -  3’~Qg”N)  log  log  N  (78) 

and  use  (76)  in  the  form 

Vi(g(N)  (79) 


From  n  we  may  then  find 

CN'  Po(V  “>  •  2 

and  round  to  the  nearest  integer.  To  obtain  more  accurate  results  n may  be  computed 
from  the  rounded  value  of  c^  and  n  could  then  be  found  from  (70)  or  (71). 

Table  2  shows  that  (76)  leads  to  good  results  for  the  two  previously  discussed  typical 
examples . 

As  a  general  conclusion  of  the  many  numerical  comparisons  which  have  been  carried  out 
we  may  state  that  the  asymptotic  formulas  give  sufficiently  good  approximations  to 
the  optimum  sampling  plans  for  most  practical  purposes.  If  one  wants  to  be  sure  to 
find  the  optimua  plan  one  may  start  from  the  approximation  and  compare  the  costs  of 
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this  plan  with  the  costs  of  suitably  chosen  neighbouring  plans  thus  finding  the 
optimum  one  by  trial  and  error. 


Table  2. 

Comparisons  of  exact  and  approximate  sampling  plans  computed  from  (76)  . 


p  =  p  =  0.010.  w  =  0.05. 
r  s  ’  2 


it 

O. 

0.006, 

p2~ 

0.040. 

pl= 

0.006, 

P2= 

0.020. 

Approx . 

Exact 

Approx. 

Exact 

N 

"n 

CN 

n 

c 

°N 

CN 

n 

c 

300 

Accept 

50 

1 

500 

20 

0 

60 

1 

i 

700 

55 

1 

65 

1 

1000 

90 

2 

115 

2 

1 

2000 

165 

3 

170 

3 

3000 

210 

4 

220 

4 

Accept 

Accept 

5000 

265 

5 

275 

5 

180 

3 

250 

4 

i 

7000 

300 

5 

285 

5 

320 

5 

335 

5  ! 

10000 

345 

6 

335 

6 

465 

7 

505 

7 

20000 

420 

8 

395 

7 

755 

10 

760 

10 

30000 

470 

0 

450 

8 

930 

12 

930 

12 

50000 

525 

9 

505 

9 

1145 

14 

1105 

14 

70000 

565 

10 

560 

10 

1290 

16 

1275 

16 

100000 

610 

11 

610 

11 

1445 

18 

1445 

18 

200000 

690 

12 

670 

12 

1750 

22 

1705 

21 

The  formulas  (72)  and  (77)  have,  however,  the  serious  drawback  from  the  point  of 
view  of  application  that  the  constants  5  and  d  are  rather  hard  to  compute.  The 
asymptotic  formulas  have  therefore  in  the  following  only  been  used  to  derive 
relationships  between  sampling  plans  under  variation  of  the  parameters.  It  i3  to 
be  expected  that  these  relationships  will  prove  to  be  rather  accurate  in  view 
of  the  good  approximation  demonstrated  above. 

According  to  (62)  we  have  for  the  optimum  plans  that  the  average  decision  loss 
asymptotically  is  constant,  i.e.  R  -  n  ~  1/c^.  For  small  N  this  gives  an  upper 
limit  to  the  decision  loss  but  the  formula  is  not  of  much  value  as  an  approximation. 

£ 


Averaoa  doci'ion  loss  as  function  c 


sions  os  functions  of  lot  s': 
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Fig.  4  sketches  for  the  two  previously  considered  examples  R  -  n  as  function  of  N. 
The  discontinuities  correspond  to  changes  in  c;  each  time  c  is  increased  by  1 
n  increases  approximately  by  p  and  R  -  n  decreases  with  the  same  quantity.  The 
asymptotic  result  corresponds  to  the  mid-points  of  the  intervals.  It  will  be  seen 
that  the  asymptote  is  nearly  being  reached  for  N  =  100,000  in  the  case  with 
P2/P1  *  6.7  but  not  for  Pj/p^  B  3.3. 

For  small  N  a  useful  upper  limit  to  the  average  decision  loss  may  be  obtained  by 
noticing  that  R  <  N?2  If  an  optimum  plan  exists  and  the  alternative  is  acceptance 
without  inspection. 

According  to  (63)  and  (64)  the  probabilities  of  wrong  decisions,  Q(p^)  and  P^)* 
are  asymptotically  inversely  proportional  to  N.  Fig.  5  sketches  for  the  two  examples 
Q(p^)  and  P(P2)  as  functions  of  N.  The  asymptotic  formula  gives  a  reasonable 
approximation  to  P(P2)  in  both  cases,  whereas  the  approximation  to  Q(p^)  is  rather 
poor,  particularly  for  the  case  P2/P1  =  3.3.  The  discontinuities  resulting  from 
changes  of  c  are  very  pronounced  for  Q(p^). 


6.  Proportional  change  of  (pf,  pg,  p2)  for  fixed  w2* 


We  shall  first  study  the  asymptotic  formulas  for  all  "quality  levels"  tending  to 
zero  with  the  same  speed.  Introducing  the  auxiliary  quantities 


P 


8 


(80) 


we  find  for  p  — >  0  and  fixed  (p  ,  p. ,  p9,  w_) 


s'  Hl' 

> 

w2(p2’1) 
ap  — >  I  In  . 

\  w^l-pj). 


(02  -  Pl)  =  v 


Pp,  — >  In  —  /  (p,  -  p. ) 


Pp< 


Po/Pi  — >  1/PQ  -  PQ, 


■Pp/Pr  — >  Pc  ln  TT  +  (pi  *  po) 


cp*, 


i  B  1  or  2; 


and 


fpr  ->  exP(<P*(a0+  I4))  }  (p  -p'ko.-p  )  °*Pf  (P0-Pt)(a0+  |°)) 

*  o  1=1  s  m  10 


“  exp  (-6q), 

Ot  1 

where  in  the  last  expression  -a  has  been  replaced  by  —  +  —  as  in  (75). 
Inserting  these  results  into  (69)  and  (72)  we  find 
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and 


ncpr  ac  +  Po(c  +  2>  *  no(c) 


'n(N  p  )  — >  cp*n  (c)  +  -  In  n  (r)  +  5 
'  crr'  o  2  o  o 


=  In  No(c). 


It  follows  that  for  small  p  wc  have  approximately 


n  -on  (c)/p 
c  o  rr 

and 

N  ~  N  (c)/p 
c  o  r 

where  nQ(c)  and  Nq(c)  are  independent  of  p^_,  i.e.  n  and  N  vary  inversely  proportional 
to  pr  for  Riven  c. 

Suppose  that  the  optimum  sampling  plans  have  been  tabulated  for  a  small  value  of 
p  ,  P  •“  0.01  say,  and  certain  values  of  (p  ,  p. ,  p_,  w_).  The  above  result  may  then 
be  used  to  find  the  optimum  plans  for  Xp^,  say,  from  the  plans  in  the  given  table. 
Denoting  the  quantities  required  by  nc(^Pr)  and  ^(Xp^)  we  have  for  given  c 


n  (Xp  ) 
c  v  rr  ' 

~  nc(pr>/>- 

(81) 

Nc(M>r> 

-  Nc(pr)/\, 

(82) 

i.e.  we  have  found  the  following  important 'proportional i ty  law"; 

The  optimum  sampling  plan  corresponding  to  (N,  Xp  ,  Xp  ,  Xp  ,  Xp9 ,  w  )  j_s 
approximately  equal  to  (n#/X,  c*)  where  (n*,c*)  is  the  plan  corresponding  to 
(N*,  Pr,  Ps,  Pp  P2,  w2)  with  N*  =  NX. 


The  theorem  has  been  illustrated  in  Fig.  6  which  shows  that  the  approximation  holds 
good  also  for  quite  large  values  of  p^. 

This  theorem  greatly  enlarges  the  field  of  application  of  the  two  master  tables. 

The  table  with  p^  =  0.01  may  be  used  for  X  <  5  and  the  table  with  p  -•  0.10  for 
0.5  <  X  <  2,  in  that  way  covering  ail  cases  with  pf  <  0.20  which  :s  the  domain  of 
practical  interest. 

A  large  number  of  numerical  investigations  has  shown  that  the  proportionality  lav; 
gives  rather  accurate  results.  The  value  of  c  found  will  seldom  deviate  more  than 
1  from  the  correct  value.  For  X  >  1  the  formula  will  normally  tend  to  give 
too  large  a  value  of  c  and  for  X  <  1  too  small  a  value. 


r 


/ 

\ 

* 


I 


Relation  between  lot  size  and  acceptance  number  by  proportional  change  0f  Pr'  Ps*  P1»  ^or  ^!Xe<i 
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Table  3. 

Comparisons  of  exact:  samplins  plans  for  p  ~  p  =  0.030,  p.  "  0.018 
p2  B  0.060,  w 2  *  0.05,  and  approximate  plans  derived  from  the  master 
tables  by  the  proportionality  law. 


N 

Exact 

n  c 

Derived  from 
pr  -  0.01  (X  =  3) 

N*-3N  n*/3 

c* 

Derived  from 

p  -  0.10  (X  »  0.3) 
r 

N*^0.3N  nV0.3  c* 

Accept 

3000 

Accept 

300 

Accept 

110 

5 

6000 

110 

5 

600 

115 

5 

3000 

140 

6 

9000 

165 

7 

900 

145 

6 

225 

9 

15000 

225 

9 

1500 

200 

8 

7000  | 

255 

10 

21000 

255 

10 

2100 

255 

10 

■ 

310 

12 

30000 

310 

12 

3000 

285 

11 

1 

395 

15 

60000 

395 

15 

6000 

365 

14 

1 

455 

17 

90000 

455 

17 

9000 

425 

16 

1 

510 

19 

150000 

540 

20 

15000 

480 

18 

| 

570 

21 

210000 

570 

21 

21000 

535 

20 

10C000 

625 

23 

- 

- 

- 

30000 

565 

21 

200000 

710 

26 

- 

- 

- 

60000 

675 

25 

The  example  in  Table  3  shows  the  derivation  of  sampling  plans  with  a  break-even 
quality  of  p^;  -  0103,  partly  from  the  first  master  table  using  X  *  3  and  partly  from 
the  second  using  X  =  0. 3.  Both  results  are  remarkably  close  to  the  exact  solution, 
see  also  Fig.  6. 


Consider  now  the  inverse  formula 
From 


we  find 


d  +  log 

nN(Xpr) 


(79). 

p  — >  0  log  e  -•  —  log  cp*  «  d 
rr  o  2  L  o 

VV  (  _  _1 - j  „  log  ,x 

X  K 3  log  N'  Xcp(pr) 


(83) 


where  cp(pr)  denotes  the  value  of  <p  for  the  given  (basic)  set  of  parameters.  This 
formula  shows  how  the  sample  size  for  a  given  lot  size  changes  with  the  "quality 
level".  This  result  is,  however,  not  as  accurate  as  the  previous  one  for  small  N 
and  it  is  neither  as  convenient  for  use  in  connection  with  the  tables. 


An  example  has  been  given  in  the  following  table  for  N  =  50,000,  p  «=  p  =  O.OJ.D, 

r  s 

p^  *  0.006,  P2  ■  0. 040,  and  W2  ~  0.05. 

Comparisons  of  exact  and  approximate  sampling  plans  derived  from  (83). 


Exact 

Approximation 

X 

n 

c 

n  c 

0.1 

1850 

3 

2330  4 

0.3 

1300 

7 

1210  7 

1.0 

505 

9 

- 

3.0 

205 

11 

210  11 

10.0 

63 

12 

78  15 
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-L  Change  of  pg  for  fixed  (pf, px, P2 , v2) . 

The  master  cables  contain  sampling  plans  for  p  -  p  only,  because  a  simple  and 

S  IT 

rather  accurate  rule  exists  for  deriving  plans  for  p  4  P  from  c^c  tabulated  ones. 

S  IT 

From  (69)  and  (72)  it  will  be  seen  that  p  influences  N  only  through  &.  Writing 

s  c 


P  -  P  “  P  -  P  +  w.  (p  -  p.) 
s  rm  rs  rr  1  rr  1 


p  -  p 

rs  rr 


=  w  (p  -  p  )(  1  +  — : - r 

l  ^r  wl(Pr"  Pp 


it  follows  from  (72)  than 


P  -  P 

1  r  s  rr 

log  N  (p  ,p  )  =  log  N  (p  ,p  )  +  log  1  +  — - 

c  r  rs  &  c  rr'Y  ^  wl'Pr’  P1  - 


or 


say,  where 


N 


c(pr,p£ 


)  =  N(p,p)/X, 
/  cv*r7  r'  s' 


X  = 
s 


P  ~  P  s 
s  r 

wi (pr “  pi>/ 


-1 


(84) 


^85) 


We  have  thus  proved  the  following  theorem: 


The  optimum  sampling  plan  corresponding  to  (N,  p  ,p  ,  p.,p0,w  )  is  approximately  equal 
*  *  -*  r  3  1  ^  ^  * 
to  the  plan  (n  ,c  )  corresponding  to  (N  ,  pr,  p^,  p^,  w2^  N  =  . 

This  theorem  makes  it  possible  to  use  the  master  tables  also  for  p  4  P  if  only  N 

* 

is  replaced  by  N  .  The  error  in  c  by  using  this  procedure  will  seldom  be  more  than 
±  1.  An  example  has  been  given  in  Table  4  with 


V  (1 


2  -  1 


-1 


0.95(1  -  0.6) 


) 


0.275. 


Table  4. 


Comparisons  of  exact  sampling  plans  for  p  =  0.010,  p  =  0.020,  p  = 

r  s  i 


0.05  with  approximate  plans  derived  from  the  master  table. 


N 

Exact 

n  c 

Approximation 
N*=  0.275N  n* 

* 

c 

300  | 

i 

Accept 

83 

5 

0 

500 

,  5 

0 

138 

10 

0 

700 

10 

0 

193 

15 

0 

1000 

!  15 

0 

275 

15 

0 

2000  i 

i  60 

1 

550 

60 

1 

3000  j 

:  110 

2 

825 

110 

2 

5000  ; 

;  120 

2 

1380 

120 

2 

7000 

170 

3 

1930 

170 

3 

10000  j 

220 

4 

2750 

220 

4 

20000 

280 

5 

5500 

280 

5 

30000 

330 

6 

8250 

330 

6 

50000  | 

390 

7 

13800 

390 

7 

70000  J 

395 

7 

19300 

395 

7 

100000 

450 

8 

2  7500 

445 

8 

200000  | 

505 

9 

55000 

550 

10 

0.006,  P2  =  0.040 
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The  corresponding  ''inverse*'  formula  becomes 

Vpr'p8>  -  Vpr'pr>  +  £  <  1  ‘  riST"  )1<>g  "s'  (86> 

Using  thi3  result  for  N  *  50000  and  the  parameters  given  in  Tabic  4  we  find 

n  -  505  -  293  x  0.9291  x  0.5607  «  350 
as  compared  to  the  exact  solution  390. 

In  the  following  sections  we  shall  limit  ourselves  to  consider  cases  with  p  ~  p 

9  f 

since  we  may  always  begin  the  analysis  by  replacing  N  by  N  if  p  =f  p  .  The  "conversion 

s  r 

factor"  X^  depends  on  w  and  the  ratios  (p  ;p.),  i.e.  )i  is  independent  of  p~  and 

•j  3  1  S  “ 

the  general  quality  level. 


8.  Proportional  change  of  (Pr/Pj.Jp2)  and  change  of  . 

Consider  the  problem  of  finding  the  optimum  plans  for  an  arbitrary  set  of  parameter 
values  (pr, p^, P2>w2>  by  combining  the  proportionality  law  with  the  relation  between 
p^  and  w£  for  given  y^  and  using  the  tabulated  plans  in  the  master  table  for  para¬ 
meter  values  (prQ, P1q, P2o>w2o^ ’  say * 

* 

The  problem  is  to  determine  k  so  that  (Pr^P^^P2)  =  (Xpfo,  Xp^Q,  Xp2q)  anc* 

(pro>  plo,P2o,W2^  givc  the  same  value  of  7 2  as  (Pro> pi07 P2o,W2o*  *  For  thiG  valuG 
of  X  we  may  find  the  plans  for  (pf,  Pj_>  P2,w2)  from  the  plans  for  (Pro;  pio,p2o,W2^ 

by  means  of  the  proportionality  law,  and  the  plans  for  (Pro>  pi0,p2o,W^  are 

identical  to  the  plans  for  (p  ,p,  ,  p«  .w^  ).  (It  will  be  noted  that  p  /p  is 

ro7  lo7  2o7  2o  ro  ro 

identical  to  the  function  defined  by  (46)) . 


Since  the  value  of  y 2 
the  equation  72q=  ^2 


is  the  same  for  (p  , p^ p2,w2)  and  (p  , 
for  the  determination  of  X,  i.e. 


w2o(p2q-  Pro) 
Ulo(pro-  Plo> 


W2(P2~  Pr) 

wi(pr-  p!> 


lo 


’p2o’v2> 


we  have 


Introducing  p2q= 


P2/x  and  P^0=  P^/X  we  find 


(87) 


For  the  master  table  with  p  =  0.01  and  w_  =  0.05  the  result  is 

ro  lo 


k  *  100(p2+  1972P1)/0  +  19>2). 


(88) 
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For  the  other  master  table  (prQ  *  0. 10  )  th®  factor  100  should  be  replaced  by  10. 

The  results  of  sections  6-8  may  be  combined  to  the  following  theorem: 

The  optimum  sampling  plan  corresponding  _t£  (N,  pr,  ps,  pp  p2,  v2)  is  approximately 
equal  to  (n*/X,  c*)  where  (n*,  c *)  may  be  found  in  the  master  table  for  N*  =  NXgX, 
P^c  ~  p^/X,  and  p2o  =  p 2 /X,  the  conversion  factors  being  equal  to 


X 

s 


Ps~Pr  \ 
vl(pr-pl>/ 


and 


X  =  100(p2  +  19  7 2  Pl)/(1  +  19  72), 


w2 (p2"pr* 

VvV  ’ 


for  the  0.01-table,  100  being  replaced  by  10  for  the  0. 10-tablc. 


By  means  of  this  theorem  it  is  rather  easy  to  find  the  optimum  plan  corresponding 
to  an  arbitrary  set  of  parameter  values  if  only  p^/X  and  p2/X  fall  within  the  range 
of  arguments  in  the  master  tables.  If  that  is  not  the  case  the  method  given  in  the 
next  section  may  be  used. 

Usually  Pj/X  and  p2/X  will  not  be  equal  to  the  arguments  used  in  the  master  tables. 
One  might  then  interpolate  but  this  is  hardly  worth  while  since  the  arguments 
in  the  table  have  been  chosen  in  such  a  way  that  by  rounding  to  the  nearest  argument 
the  rounding  error  will  ordinarily  be  less  than  10  %. 

If  one  wants  to  be  sure  to  obtain  a  sufficiently  large  sample  the  value  of  p^/X 
should  be  rounded  up  and  the  value  of  p2/X  rounded  down. 

As  an  example  consider  the  problem  of  finding  the  sampling  plans  for  (pr,p^,p2,w2) 

»  (0.03,  0.01,  0.07,  0.08)  and  p  -  p  .  Since  p  <  0.05  say,  wc  choose  to  use  the 

S  r  IT 

the  0.01-table.  From 

*2  -  4  £f  =  °-m>  19  ^2  =  3-31’ 

we  find  X  *  (7  +  3. 31 )/  (1  +  3.31)  =  2.39,  p^X  -  0.01/2.39  *  Q.0042,and 

p0/X  r  0.07/2.39  *  0.029.  The  master  table  should  thus  be  entered  with  p,  =  0.004 
Z  rlo 

and  p2q  =  0.030.  For  N  =  2000,  say,  we  find  N*  =  4780  and  (n*,c*)  =  (210,  3)  leading 
to  (n,c)  -  (210/2.39,  3)  -  (90,3)  which  is  the  correct  solution. 

If  w2  =  0.02  instead  of  0.08  we  find  similarly  X  =  4.38,  p^/X  =  0.0023s;  0.0025  and 
P2/X  =  0.0160  =  0.0150.  For  N  -  2000  we  get  N*  =  8760  leading  to  acceptance  without 
inspection  as  the  most  economical  decision. 
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9.  Change  of  w2  for  fixed  (pf,  pg,  p1#  p^)  . 

In  the  following  we  shall  develop  a  method  for  evaluating  the  effect  of  changing 
one  of  the  five  parameters  only,  and  use  it  first  for  w2  and  then  for  Pr» 


From  (69)  and  (72)  wc  find  for  given  c 

dn 


d  a 


Ti 


og  w2  olog  w2 


ql 

1/  (w«  1  og  ■—) 
1  <12 


(09) 


and 


dlog  N 

c 

dlog  w2 


dn 


dlog  n 

+  -  - -  + 


as 


olog  w2  '  2  dlog  w2  ‘  Olog  w2 


(90) 


The  last  term  om  the  right  hand  side  is  a  rather  complicated  function  of  the 
parameters.  Tabulation  of  6  and  graphical  analysis  of  5  as  a  function  of  log  w2  has 
shown,  however,  that  at  least  for  w2  ^  0.20  and  p^  ^  0.10  (and  corresponding 
values  of  pg,  p^,  p2 )  8  is  approximately  a  linear  function  of  log  w2  with  a  slope 
depending  slightly  on  (p  ,  p  ,  p„ )  and  being  practically  independent  of  p  . 

Limiting  ourselves  to  the  case  p  =  p  we  thus  have 

S  IT 


d  5  ,  . 

diog  w2  ~  '  bi(pr  V' 

say,  where  b^fp^,  p2)  has  been  tabulated  in  the  appendix. 

Writing  8  =  &(p^,  p^,  p2,  w2)  and  putting  w2  =  0.02  and  0.20  respectively,  so  that 
so  that  Alog  w2  -  log  C., 20  -  log  0.02  -  1,  an  approximation  to  dh/dlog  w2  may  be  found 
as  S(Pr>  P^j  P2>  0.20)  -  8(pr,  Pj^,  p2,  0.02).  This  approximation  has  been  computed 
for  both  =  0.01  and  0.10,  and  finally  the  average  of  the  two  has  been  taken  as 

-br 


For  small 


p  wc  also  hove 
r 

9 


W  (p2 


Px)  =  b2(p). 


The  values  given  for  b2  have  been  computed  as  averages  of  9/ (log 

p  -  0.01  and  p  «=  0.10. 
r  r 

For  large  n  we  have  that  (log  e)/2n  is  small  as  compared  to  9  and 
disregard  the  second  term  on  the  right  hand  side  of  (90).  We  then 


ql 

—)  for 
q2 

we  shall  therefore 
have  approximately 


/ 


dlog  Nc  b2(p) 
dlog  w2  w1 


-  b^p^  p2) 
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which  gives 

W  -  AW 2'b^ 

where  A  denotes  a  constant  of  integration.  Changing  from  to  Xw2  we  get 


Nc(Xw2)  -=  Nc(w2)/f1(X) 

where 

£l<»  -  x”1  bj(l  -  <x  -  i>JJj) 


(91) 


(92) 


From  (69)  we  further  have 

nc^Xw2^  “  nc(w2>  +  81^X* 


where 


gj(X) 


(p2-p1)Pr. 


(93) 


(94) 


For  convenience  f^  and  g^  have  been  written  as  functions  of  X  only,  even  if  they  both 
depend  also  on  other  parameters.  The  function  f^  which  will  be  called  the  conversion 
factor  for  N  due  to  a  change  in  w2  has  been  tabulated  in  the  appendix  for  w2  =  0.05  as 
a  function  of  (X,p^,p2).  The  function  g^  which  gives  the  correction  to  n  due  to  a  change 
in  w2  has  similarly  been  tabulated  in  the  appendix  as  function  of  (x,p^,p2)  for  v»2=  0.05 
and  pr*  0.01.  Values  of  this  function  for  other  values  of  p^  may  be  obtained  as  g^/100pr 
where  g^  represents  the  tabulated  values. 

The  above  results  may  be  formulated  as  the  following  theorem: 


The  optimum  sampling  plan  corresponding  to  (N,p  ,p  , p. , p_ ,  \w_) ,  p  *  p 
equal  to  (n*+  g^(>),c*)  where  (n*,c*)  is  the  plan  corresponding  to  (N* 


,  is  approximately 
>pr’Ps,IVP2'W2* 


with  N*=*  Nf1(\) . 


The  theorem  has  been  illustrated  in  Fig.  7. 


This  theorem  enlarges  the  field  of  application  of  the  two  master  tables  with  respect  to 
values  of  w2  in  a  similar  manner  as  the  law  of  proportionality  does  with  respect  to  the 
other  parameters.  The  results  of  using  the  approximation  have  been  compared  with  the 
exact  solutions  in  a  large  number  of  cases  and  the  deviations  found  between  the 
approximate  and  the  correct  value  of  c  have  never  exceeded  1  for  X  <  4.  There  is  a 
tendency  for  the  approximation  to  give  too  small  a  value  of  c  for  X  >  1  and  too  large 
a  value  for  X  <  1,  in  particular  for  small  N. 


Relation  between  lot  size  and  acceptance  number  by  change  of  W2  for  fixed  (Pf, 
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It  should  be  noted  that  the  formula  breaks  down  in  some  cases  for  small  N.  Let 

denote  the  largest  N  for  which  acceptance  without  inspection  is  cheaper  than 

sampling  inspection  for  the  master  table  used.  If  \  >  1  and  Nf,(X)  *  N*  <  N  then 

1  a 

the  formula  does  not  lead  to  a  sampling  plan  even  if  there  may  exist  a  plan  which 
for  \w^  is  cheaper  than  acceptance  without  inspection.  Similarly,  for  X  <  1  and 
Nf^(X)  “  N*  >  there  may  be  3ome  cases  where  the  approximation  formula  leads  to 
a  sampling  plan  even  if  the  cheapest  solution  is  acceptance  without  inspection. 


/in  example  has  been  shown  in  Table  5.  The  approximation  is  remarkably  good.  Since 
N  *  74  the  approximation  formula  leads  to  acceptance  without  inspection  for  all 

a 

N  <  57.  Sampling  plans  cheaper  than  acceptance  without  inspection  do,  however, 
exist  for  12  £  N  $  57. 

Using  the  method  of  section  8  we  find  y^~  0.833,  X  =  0.804,  p^/X  =  0.0075,  and 
P2/X  e  0.050.  Since  p^/X  falls  outside  the  range  of  arguments  in  the  master  table 
the  method  does  not  apply.  Using  p^/X  *  0.007  gives,  however,  a  rather  good  approxima¬ 
tion. 


Table  5. 


Comparisons  of  exact  sampling  plans  for  p  =  p  =  0.010,  p.=  0.006, 
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P2=  0.040,  w^*  0.10,  and  approximate  plans  derived  from  the  toaster 
table,  f^  1.29,  gj=  20. 


N 

Exact 

n  c 

N*=  1.29N 

Approxima 

n*+  20  c 

50 

15 

0 

65 

Accept 

70 

20 

0 

90 

25 

0 

100 

25 

0 

129 

30 

0 

200 

35 

0 

258 

35 

0 

300 

75 

1 

387 

75 

1 

500 

85 

1 

645 

85 

1 

700 

130 

2 

903 

130 

2 

1000 

140 

2  | 

1290 

140 

2 

2000 

240 

4 

2500 

240 

4 

3000 

250 

4  i 

3870 

250 

4 

5000 

305 

5 

6450 

300 

5 

7000  1 

355 

6 

9030 

355 

6 

10000  ; 

405 

7 

12900 

405 

7 

20000 

465 

8 

25000 

465 

8 

30000 

520 

9 

38700 

520 

9 

50000 

575 

IO 

64500 

575 

10 

70000 

630 

11 

90300 

630 

11 

100000 

635 

11 

129000 

635 

11 

200000 

740 

13 

- 

- 

- 
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From  the  inverse  formula  (79)  we  get 


and  consequently 


or 


,nN 


d  log  w„ 


b. (Pi) Po)  1 

1  1  ~  (l  -  - - -  ) 

\  1  -1  ^ _ M  / 


(P 


3  log  N 


Vw2>  “  A  +  -  d  -  jISTn)  log  w2 


V'V  =  W  +  ~~  U  '  3  log"M>  108  *  • 


(95) 


This  shows  that  the  difference  between  n^w^)  and  f°r  gi-ven  N  is  proportional 

to  log  X.  This  formula  is,  however,  not  as  accurate  as  (93)  for  small  N. 

An  example  has  been  given  in  the  following  table  for  N  =  20,000,  p  =  p  -  0.010, 
p^  0.006,  p^=  0.020,  and  w^3  0.05,  which  gives  b^=  0.59  and  1/cp  *  1100. 

Comparisons  of  exact  and  approximate 
sampling  plans  derived  from  (95)  . 


Exact 


Approx. 


and 


100  w^  X 

n 

c 

n 

c 

2.5  0.5 

540 

8 

580 

9 

5.0  1.0 

760 

10 

- 

- 

10.0  2.0 

980 

12 

940 

11 

20.0  4.0 

1130 

13 

1120 

13 

Pr=  Pr  for  fixed  (p1?p^,w0) 

• 

(72)  we  find  for  given  c 

da 

/ 

1 

+ _ L 

_\ 

* 

1 

1 

H-* 

d  log  Pr  '  '3  log  pr 

’Pr 

Pr-  PX 

p2‘ 

V, 

/  '  q2/ 

>1  log  N  dn 

c  c 

1  d  log 

n 

_c  . 

do 

d  log  Pr  ~  r  d  log  pr 

2  d  log 

■  T 

Pr 

5” 

lo8  Pr 

(96) 


(97) 


Numerical  investigations  show  that  -  for  w^  <  0.20,  p^  <  0.20,  and  p^=  p  -  6  is 
approximately  a  linear  function  of  log  p^  with  a  slope  depending  on  u  and  being 
practically  independent  of  "the  level  of  (p^,P2)"  and  of  if  only  p^  does  not 
come  too  close  to  p^  or  p£,  i.e. 


i>6  ....  1/5  _  -1/5 

ISTT  ^  b3(fl)  £or  plp  <  pr<  P2C 
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say,  where  b^(p)  has  been  tabulated.  (Another  limitation  of  no  practical  importance 
ia  that  p2  must  not  be  too  close  to  1).  An  approximation  to  d8/dlog  pf  nay  be 

computed  as  the  corresponding  difference  -  quotient  setting  p  -  ;.p  and 

-  1/5  r  1 

Pr  “  P2P  respectively.  This  has  been  done  for  *=  0.05  and  for  the  "ntandard" 

values  of  p^  and  p2,  partly  at  the  1  7,  and  partly  at  the  10  %  level.  The  value  of  b^ 

given  in  the  table  is  the  average  of  the  two  values  found. 

Proceeding  as  in  section  9  we  have  approximately 


dlog  N 

c 

Pr 


b,  (p)p  ( - + 

2  ^/rr  p  - p, 

r  1 


P2-Pr 


)  + 


b3(p) 


which  on  integration  gives 


Nc(Pr) 


and 


N  (Xp  ) 
c  r 

where 

f2  (X) 

From  (69)  we  further  have 


Nc(Pr)/f2(X) 

,"b3^  (Oj-1)(X-P1)^b2 

x  Vo-ppCPj-x)/' 


nc  (>-pr)  ”  nc<Pr)  +  x2  (>-) 


where 


g2(x) 


at  In 


(i-p1)(p2-x) 

(p2-l)(X-p1) 


(P2  '  Pi>Pr- 


(98) 

(99) 


(100) 


(101) 


The  conversion  factor  for  N  due  to  a  change  in  pf,  f2(X),  has  been  tabulated  in  the 
appendix  as  function  of  (X,  pj ,  p^),  and  the  correction  to  n  due  to  a  change  in 
Prt  82(X),  has  been  tabulated  as  function  of  (X,  p^  p2)  for  pr  =  0.01.  Values  of 
g2  (X)  for  other  values  of  p^  may  be  found  from  the  tabulated  ones  by  dividing  by  lOOp  . 


The  above  results  may  be  formulated  as  the  following  theorem  ; 


The  optimum  sampling  plan  corresponding  to  (N,  Xp  ,  Xp  , 

I*  s 

approximately  equal  to  (n*  H-  g2(X),  c*)  where  (n*,c*)  i£ 
(N*,  Pr,  pg,  PL,  P2,  w2)  with  V*  =  Nf 2  (X) . 


pl'  P2'  W2^  pr  =  V 
the  plan  corresponding  to 


With  the  given  set  of  tables  this  theorem  is,  however,  not  as  important  in  practice 
as  the  previous  ones,  because  the  tables  contain  the  optimum  plans  for  so  many 
combinations  of  (Pr>  P^>  P2)  that  an  adjustment  of  the  relative  position  of  p  vi thin 
the  interval  (p^,  p2)  will  seldom  be  felt  necessary  from  a  practical  point  of  view. 


In  table  6  an  example  has  been  shown  of  the  effect  of  changicg  p  ■  p  from  0.010 

r  8 

to  0.020  within  the  interval  (p^,  p2)  ■  (0.006,  0.040). 


Table  6. 


Comparisons  of  exact  sampling  plans  for  p  ■  p  ■=  0.020,  p.  “  0.006, 

IT  S  A 

P2  “  0.040,  w2  ■  0.05  with  approximate  plans  derived  from  the  master 
table.  X  -  2,  f2  «  0.52,  g2  «  -  55. 


N 

Exact 

n  c 

.  Approximate 

N*  -  0.  52N  n*-55 

c* 

2000 

Accept 

1040 

60 

2 

3000 

75 

2 

1560 

110 

3 

5000 

130 

3 

2600 

165 

4 

7000 

180 

4 

3640 

170 

4 

10000 

230 

5 

5200 

225 

5 

20000 

290 

6 

10400 

280 

6 

30000 

345 

7 

15600 

335 

7 

50000 

400 

8 

26000 

390 

8 

70000 

450 

9 

36400 

445 

9 

100000 

460 

9 

52000 

450 

9 

200000 

565 

11 

104000 

555 

11 

From  the  inverse  formula  (79)  we  get 

S  .  .  „  .  _i _ , 

Slog  pr  <p  '  3  log  N  ' 

which  leads  to 

VXpr>  -  n„(Pr)  -fd-  TTSTn  >  108  X-  (102  > 

An  example  of  the  application  of  this  formula  has  been  given  in  the  following  table 
for  N  =  50,000,  Pr  “  Pg  *  0.010,  X  **  0. 5  and  2.5,  Pj^  -  0.002,  p2  ■  0.040,  and  w2  «  0.05, 
which  give  -  1.09  and  l/<p  «  177. 

Comparisons  of  exact  and  approximate  sampling  plans  derived  from  (102) 

Exact  Approximate 


pr 

X 

n 

c 

n 

c 

0.005 

0.5 

350 

4 

370 

4 

0.010 

1.0 

315 

4 

- 

- 

0.025 

2.5 

250 

4 

245 

4 

11.  Change  of  all  parameters. 

I 

The  results  of  the  preceding  sections  may  be  combined  intc  a  •'chain  formula"  of  the  type 
Nc(Xpr,  PsXpr,  \9V  Xp2,  XjW2)  *  Nc<Pr>  Pr>  Pp  P2>  w2)/XsfiX  (103) 

and 

nc(Xpr,  PsXpr,  XpL,  Xp2,  Xxw2)  -  (nc(pr,  pr,  p^  p2,  w2)+  g^/X 


(104) 


where 


-1 


-  42  - 

,  i,  ,  p.~‘ 

»  '  V  P-x^Ki-o,) 

and  g1(X1)  being  defined  by  (92)  and  (94)  for  =  p1/pr  and  p2  =  P?/Pr' 

In  the  master  tables  p  ■  p  =0.01  (or  0.10)  and  w~  -  0.05  have  been  used  as 

r  8  4 

reference  values.  What  has  been  denoted  by  X  and  X^  in  the  above  formulas  become  lCGp^ 
(or  lOp^)  and  20w2  if  pf  and  w2  denote  the  values  for  which  the  optimum  plan  is 
required. 

We  thus  get  the  following  rule  for  using  the  master  table  with  p^  =  O. 01 : 

The  optimum  plan  for  (N,  p  ,  p  ,  p.  ,  p„ ,  w„)  with  p  <  0.05  is  approximately  equal  to 

r  S  1  /  t  IT 

((n*  +  g^)/100pr,  c*)  where  (n*,  c*)  may  be  found  by  entering  the  master  table  with 

N*  -  H(100pr)f1(20w2,  Pl>  p2j (l  +  )  ,  (105) 

P1  “  pi/pr^  p2  "  p2^pr'  SM  "  gj^OWj,  pl>  p2^  the  arguments  for  (px,  p2) 
in  the  masxer  .table  being  (p^/lOO,  p2/100). 

For  0.05  <  pr  <  0.20  the  master  table  with  p^  -  0.10  should  be  used  accordingly. 

If  (p^/lOO,  p2/100)  are  not  to  be  fouiad  in  the  table  then  use  the  "nearest"  argument 
or  interpolate.  One  may  also  use  the  results  in  section  10  to  change  pf  in  the  master 

table  so  that  the  relations  between  (p^,  p^,  p2)  in  the  table  become  closer  to  the 

ones  for  which  the  sampling  plan  is  required.  From  a  practical  point  of  view,  however, 
the  master  tables  combined  with  the  rule  above  will  normally  suffice. 

An  example  has  been  given  in  Table  7.  The  conversion  factor  for  N  is  found  as 

3f  L  (2,  0.6,.  4.0)  /  (1  +  o79o-°-x-  jj)  -  3  *  1.29/3.78  =  1.02. 


The  agreement  between  the  approximate  and  the  exact  solution  is  very  good. 

Using  instead  the  method  of  section  8  we  get  X  =  1/3.78  =  0.265,  X  =  2.40, 

s 

p^/X  *=  0.0075,  and  P2/X  =  0.050,  i.e.  N*  =  0.636N  and  n  =  n*/2.40.  Since  the  master 
table  does  not  contain  the  argument  0.0075  wc  may  as  an  approximation  use  0.0070 
which,  however,  will  tend  to  give  too  small  values  of  c. 

The  corresponding  inverse  formula  takes  the  form 


WWW  "  To0^[no  +  l  {l°8«°°pr>  +  Vog(20w2) 

-  lo/l  +  — ;-S  -rr  yi  (1  -  y-p — -) 

\  w1(pr-p1)  J)  3  log  N  J 


(106) 
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where  nQ  denotes  the  sample  size  to  be  found  in  the  master  table  for  =  0.01, 
pf  *=  P]/Pr>  p2  =  P2^Pr'  corresponding  to  the  given  lot  size  N. 


Table  7. 

Comparisons  of  exact  sampling  plans  for  p  *  0.030,  p  *  0.060, 
p^  =  0.018,  P2  *  0.120,  =  0.10  and  approximate  plans  derived  from 

the  master  table  for  Pr  “  Ps  =  0.010,  p^  =  0.006,  P2  “  0.040,  W2  =  0.05  • 


N 

n 

Exact 

c 

N*  =  1.02N 

Approximate 
(n*+  20)/3 

c* 

50 

Accept 

51 

Accept 

70 

5 

0 

71 

Accept 

100 

5 

0 

102 

10 

0 

200 

10 

0 

204 

10 

0 

300 

10 

0 

306 

25 

1 

500 

25 

1 

510 

25 

1 

700 

30 

1 

714 

30 

1 

1000 

45 

2 

1020 

45 

2 

2000 

65 

3 

2040 

65 

3 

3000 

80 

4 

3060 

80 

4 

5000 

100 

5 

5100 

100 

5 

7000 

100 

5 

7140 

100 

5 

10000 

120 

6 

10200 

120 

6 

20000 

135 

7 

20400 

155 

8 

30000 

155 

8 

30600 

155 

8 

50000 

175 

9 

51000 

175 

9 

70000 

190 

10 

71400 

195 

10 

100000 

210 

11 

102000 

210 

11 

200000 

225 

12 

204000 

230 

12 

As  an  example  consider  the  determination  of  n  for  N  =*  50,000  and  the  parameters  given 
in  Table  7.  The  value  of  1/ cp  is  293  for  p^  =  0.006  and  p2  =  0.040,  and  (1  -  3  log  tP 
*  0.929,  so  that  we  find 


n  =  3<505  +  293 (log  3  +0.61  log  2  -  log  3.78)0.929) 

=  3(505  +  23)  =  176 

in  agreement  with  the  (rounded)  exact  solution,  n  =  175,  given  in  Table  7. 


12,  Efficiency. 

In  a  previous  paper  {_  6  _]  it  has  been  proposed  to  define  the  efficiency  of  a 
sampling  plan  as 

e(N,  n,  c)  -  Rq(N)/R(N,  n,  c)  (107) 

where  R(  ,(*)  denotes  the  costs  of  the  optimum  plan  and  R(N,  n,  c)  denotes  the  costs 
of  the  plan  in  question. 
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We  shall  first  discuss  the  efficiency  of  a  sampling  plan  on  the  asstu^tion  that 
the  aptimun  relationship  between  n  and  c  has  been  used  so  that  the  loss  in 
efficiency  is  due  to  using  a  wrong  relationship  between  N  and  n.  Looking  at  Fig.  2 
it  will  be  seen  that  it  does  not  matter  much  whether  we  use  ♦’ha  value  of  c  giving 
the  absolute  minimum  of  R  or  a  neighbouring  value  of  c  provided  n  is  chosen  such 
that  a  (relative)  minimum  of  R  is  obtained. 


For  a  given  set  of  parameters  let  (n0>cQ)  be  optimum  for  No  and  (n^,c^)  be  optimum 
for  N^.  From  (26)  it  follows  that 


R(N,  n,  c) 


n  +  (N  -  n)h(n,  c) 


where 

h(n,  c)  -  71Q(P1)  +  72p(P2>* 

Using  the  plan  (n^,  c^)  for  lot  size  Nq  (instead  of  N^)  we  find 


R<No#ni^ci)  °  +  <N0  “  n1)h(n1,  cx) 


-  nL  +  (Nq  -  n1)(Ro(N1)  -  n^/O^  -  n^.  (108) 

It  is  therefore  rather  simple  by  means  of  the  function  Rq(N)  to  evaluate  the 
efficiency  of  plans  contained  in  the  master  table  in  case  such  plans  are  used  for 
the  wrong  value  of  N.  The  resulting  efficiency  is 


R0<N0> 

e(Nb,n1^c1)  “  n1  +  (NQ-n1)(Ro(N1)-n1)/(N1  -  n^ 

Since  Rq(N)  n  +  l/cp^  we  have  asymptotically 


e(No,Vci)  ^ 


Vni 

Nrni 


(109) 


(l.’O) 


Introducing  n^  *■  (In  N^)/^  +  o(ln  Nq)  and  considering  as  an  arbitrary  function 
of  Nq,  n^  =  8(NQ)/<P0  say,  we  find 

e(N,  n^Cj)  ^  (In  N)/(g(N)  +  Ne~s(N))  (111) 

for  n^  -  o(N),  which  is  the  result  given  without  proof  in  [  6  J. 

For  g(N)  =  X  In  N  we  get  e  — >  1/X  for  X  ^  1  but  e  — ■>  0  for  0  <  X  <  1,  i.c.  if  we 

use  a  semilogaritmic  relationship  between  n  and  N  differing  from  the  correct  one 

then  it  is  important  to  use  too  large  a  sample. 


Ub  - 


For  g(N)  »  N^,  X  >  0,  we  get  e  — >  0.  A  more  accurate  expression  than  (111)  may 
be  found  by  using  all  three  term  of  (61)  which  leads  to 


e(N0,ni>ci)  ~ 


This  formula  is,  however,  not  of  direct  value  because  it  contains  which  is 
unknown  in  practice.  A  simple  and  practically  useful  approximation  is  the  following 


e(No,ni,ci)  ^ 


(n-n.  )' 
o  o  i 


(112) 


This  formula  will,  however,  give  too  large  efficiencies  for  small  values  of  n 
because  the  decision  loss  has  been  overestimated. 


In  connection  with  the  various  approximations  developed  in  the  preceding  sections 
is  has  repeatedly  been  stated  that  the  value  of  c  found  by  using  the  approximations 
will  normally  not  deviate  more  than  1  from  the  correct  value  (for  N  <  200,000). 

It  is  therefore  of  importance  to  know  the  efficiency  of  a  plan  for  which  lci"c0l  =  1. 

If  lci*c0l  *  a  (constant)  then  |n^-no|  =  ap  and  e  — >  1  for  Nq  — >  «>,  Expanding  the 

denominator  of  (112)  we  find  for  small  values  of  9a0 

o 

c<No'VCl>  '  <"o  +  7*)  '  <“o  +  r  +  2  <U3> 

Yo  ^o 

which  converges  rather  fast  to  1  for  nQ  — >  »  and  a  =  1.  By  means  of  the  results 
of  section  6  it  will  be  seen  that  this  asymptotic  efficiency  (as  a  function  of  cq) 
is  independent  of  the  "quality  level". 

An  example  has  been  given  in  Table  8.  The  costs  for  each  optimum  plan  have  been 
compared  with  the  costs  of  using  a  neighbouring  plan,  i.e.  c^  *  cq  ±  1.  The  efficiency 
has  been  ccspared  vitb  the  asymptotic  efficiency  found  from  (112).  It  will  be 
seen  that  the  efficiency  is  larger  than  0.90  for  c  5  2  and  that  the  asymptotic 
formula  gives  too  high  an  efficiency  for  small  N.  (N  has  been  chosen  as  the 
geometric  mean  of  the  smallest  and  largest  N  for  each  c).  This  conclusion  is 
typical  for  the  cases  investigated. 
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Table  8. 


Investigation  of  efficiency  for  sampling  plans  with  an  acceptance  number 
deviating  1  from  the  optimum*  (e*  ■  asymptotic  efficiency). 

Pr  *  pfl  ■  0.010,  Pj  »  0.006,  p2  ■  0.040,  v2  -  0.05. 


N 

7 - 

1  ni  •  Cj 

lOOe  lOOe* 

145 

10  0  53 

60  1 

72  94 

447 

60  1  126 

10  0 

115  2 

84  94 

86  95 

1010 

115  2  200 

i 

_ 

60  1  90  95 

170  3  92  96 

1900 

:'DTL  E9K9  Wii . ■ 

3350 

• 

.  ,,,  '170  3  95  96 

225  4  326  !  280  5  96  97 

1  i 

5700 

280  5  386 

225  4  96  97 

335  6  ;  97  98 

96  97 

97  98 

15800 

1 

390  7  502 

335  6 

445  8 

97  97 

98  98 

23800 

445  8  559 

390  7 

500  9 

97  98 

98  98 

42100 

500  9  616 

445  8 

555  10 

97  98 

98  98 

68300 

555  10  672 

500  9 

615  11 

98  98 

98  98 

The  conversion  formulas  and  tables  show  how  sensitive  the  solution  is  to  changes  of 

the  parameters.  A  change  of  w2  from  0.05  to  0.10,  say,  means,  that  N  has  to  be 

multiplied  by  a  factor  of  about  1.3  and  the  corresponding  n  should  be  increased  by 

about  30.  (In  most  systems  of  sampling  plans  in  practical  use  to-day  the  same 

plan  is  used  for  a  rather  large  N  interval,  the  ratio  between  endpoints  usually 

being  1.5  or  larger).  As  an  example  consider  the  case  wilh  p  =  p  *  0.010, 

r  s 

Pj  *  0.006  and  p2  *  0.040  as  shown  in  the  following  table. 
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Optimum  sampling  plans. 


V  = 
2 

0.05 

n 

CM 

* 

0.10 

N 

n 

c 

n 

c 

500 

60 

1 

85 

l 

1000 

115 

2 

140 

2 

5000 

275 

5 

305 

5 

10000 

335 

6 

405 

7 

50000 

505 

9 

575 

10 

100000 

610 

11 

635 

11 

For  most  lot  sizes  we  find  the  same  value  of  c  and  a  difference  in  n  of  about  25, 
in  other  cases  the  difference  in  c  is  1  and  the  difference  in  n  correspondingly 
larger.  It  is  immediately  clear  that  using  the  plans  corresponding  to  w2=  0.05  if  the 
true  value  of  w2  is  0.10  does  not  lead  to  an  essential  loos  in  efficiency. 

The  conclusion  is  that  even  if  the  value  of  w2  used  deviates  from  the  true  value 
by  a  factor  of  2  the  method  will  nevertheless  lead  to  a  sampling  plan  of  very 
high  efficiency. 

Similar  conclusions  may  be  drawn  for  the  other  parameters  by  studying  the  conversion 
formulas . 

The  main  reason  why  changes  of  p  and  w^  does  not  affect  the  optimum  solution  serious¬ 
ly  is  that  pQ  and  .p^  are  independent  of  pf  and  W2 . 

Since  the  most  important  relation  in  the  system  is 

c  +  j  =  PQ(n  -  ot) 


it  is  of  importance  to  know  how  pQ  depends  on  p^  and  p2 
From 


6  In  p 
_ _o 

<5  In  p^ 


V  P1 


qi1"  ~ 


>  0  and 


ct  In  p 

_ _c 

o  In  p„ 


P2~  P1 


q„ln  — 

2  ’2 


>  0 


it  follows  that  pQ  is  an  increasing  function  of  as  well  p^  as  p2  .  Furthermore  we 
have  approximately 

0  In  p  0  In  p 

- 2.  +  - 2  s  1 

d  In  Pl  In  p2 


Within  the  domain  of  variation  tabulated  the  first  term  is  on  the  average  0.35  and 
the  second  0.65. 

The  coefficient  p^a  varies  rather  slowly  with  (p^,p2). 


48 


CT\ 


Q) 

H 

JO 

H 
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It  follows  that  p  Is  known  with  a  relative  error  of  about  the  same  size  as  the 

o 

relative  errors  of  and 

If  the  choice  of  p^  and  is  doubtful  then  p^  should  be  chosen  too  large  and  p£ 
too  small  (by  about  half  of  the  percentage  error  in  p^)  because  the  two  errors 
will  tend  to  counterbalance  one  another  and  thus  give  the  correct  pQ.  The  reason 
for  bringing  the  two  parameters  closer  together  in  case  of  doubt  lies  also  in  the 
fact  that  qpo  is  a  decreasing  function  of  p^  and  an  Increasing  function  of  P2.  Since 
n  -v  (log  N)/<P0  the  proposed  rule  will  lead  to  a  larger  sample  size  than  the 
optimum  one  which  normally  gives  a  better  efficiency  than  too  small  a  sample.  However, 
with  a  wrong  pQ  (and  9^)  the  efficiency  will  tend  to  zero  because  the  average  decision 
loss  becomes  0(1^),  0  <  e  <1,  instead  of  0(1)  as  for  the  optimum  plan. 

Table  9  shows  the  efficiency  of  using  a  plan  obtained  by  entering  the  master  table 
by  a  wrong  value  of  P^>p2  or  both.  It  is  assumed  that  the  true  values  of  (p^,p2) 
are  (0.006,0.040)  and  optimum  plans  have  been  substituted  by  plans  obtained  by  using 
neighbouring  values  of  (p^,P2)  in  the  tables,  i.e.  the  relative  error  of  p^  is  177. 
and  the  relative  error  of  P2  is  12.57.  downwards  and  25%  upwards.  The  table  shows 
that  the  efficiency  in  all  cases  is  larger  than  907.  for  N  <  10,000.  For  N  =  200,000, 
however,  the  efficiency  falls  to  58%  in  the  worst  case,  i.e.  the  case  where  p^  is 
chosen  25%  too  large. 

The  results  in  the  table  support  the  statement  above  that  in  case  of  doubc  it  is 
Important  to  use  a  large  value  of  p^  and  a  small  value  of  P2> 


13.  An  example. 

Consider  now  an  example  starting  from  the  original  cost  functions.  To  show  the 
various  aspects  of  the  method  the  example  will  be  worked  out  in  more  detail  than 
is  necessary  for  routine  applications. 

Let  the  three  cost  functions  be  k  (p)  ■  23  +  35p>  k  (p)  «  16  +  35p,  k  (p)  ~  720p, 

9  IT  Q 

the  coefficients  denoting  costs  per  item  in  cents,  say,  i.e.  the  costs  of  sampling 
and  testing  is  23  cents  per  item  in  the  sample  and  the  costs  of  accepting  a  defective 
item  is  720  cents  etc.,  see  section  2. 


Let  us  further  assume  that  lots  are  generated  with  probability  w^=  0.93  from  a 
binomial ly  controlled  process  with  p^=  0.009  and  with  probability  W2=  0.07  from 
a  process  with  P2*  0.080. 


The  costs  may  then  be  described  as  in  the  following  table: 


w 

0.93 

0.07 

Average 


P  k8<P> 

0.009  23.315 
0.080  25.800 
0.014  23.489 


kr(p) 


ka(p)  kn(P>  lkr(p)"ka(p)  | 


16.315  6.480 
18.800  57.600 
16.489  10.058 


6.480  9.835 
18.800  38.800 
7.342  11.863 
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From  (12)  we  find 


pr*  (16-0)7(720-35)  =  0.0234, 


from  (28) 


and  from  (22) 


P  «  0.93  x  0.009  +  0.07  x  0.0234  =  0.0100, 

d 


p  -  (23-0) /(720-35)  =  0.0336. 
s 


To  find  the  optimum  plan  for  N  -  500  from  the  master  table  with  p  =  p  -  0.010  we 

X*  s 

first  have  to  find  the  conversion  factor  \  which  corrects  for  the  difference  between 

s 

p  and  p  ,  i . e . 

*s  rv* 


1  + 


336-234 

0.93(234-90) 


1.76. 


To  use  the  method  of  section  C  we  find  y  -  0.296,  \  -  1.97,  p^/\  =  0.0046  &  0.005, 
P2/\  “  0.041  s  0.040,  and  N*  -  1.97  N/1.76  =  1.12  N  =>  560.  From  the  master  table 
we  read  (n*,c*)  «  (60,1)  which  gives  n  «  60/1.97  -  30  as  the  optimum  sample  size. 

To  illustrate  the  method  of  section  11  we  have  to  find  the  conversion  factor  f^ 
and  the  correction  corresponding  to  the  change  from  w2“  0.05  to  0.07.  Since 

Pj*  90/234  *  0.385  sc  0.40  and  800/234  *  3.42  s  3 . 50  we  have  f^=  1.13  and  15. 

We  Chen  enter  the  master  table  with 


N*  «  N  x  2.34  x  1.13/1.76  -  1.50N  -750 

and  find  (n*,c«)  *  (60,1)  which  finally  gives  (n,c)  =  (30,1)  since  (60  +  15)/2.34  =  32. 


To  find  the  corresponding  value  of  R  we  first  compute 

0.93 (2 34- 90)/ (33 6- 100)  =  0.567 

and 

72=  0.07(800-234)/ (336- 100)  «  0.168 

which  lead  to 

R  -  n  +  (N-n)(0.567Q(p1)  +  0.168P(p2)). 

From  a  table  of  the  binomial  distribution  one  finds  for  (n,c)  *  (30,1)  that 
Q(P1)  ■  0.02982  and  P(p2)  ~  0.29579  and  consequently 


R  =  30  +  470  x  0.0666  =*  61.3. 
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The  costs  of  sampling  inspection  and  the  average  decision  losses  per  lot  are  thus 
of  nearly  the  same  size. 

Returning  to  the  original  monetary  unit  we  find 

k  -  k  =  R(k  -  k  )/500  =  1.98 
m  s  m 

and  finally 


k  «  7.34  +  1.98  =  9.32. 


We  thus  have  the  following  conclusion: 

The  quality  of  submitted  Iocs  is  such  that  on  the  average  costs  per  item  will  be 
7.34  cents  if  all  lots  are  classified  correctly,  i.c.  all  lots  from  process  No.  1 
are  accepted  and  all  lots  from  process  No.  2  arc  rejected.  To  decide  whether  to 
accept  or  reject  we  inspect  a  sample  of  30  items  at  the  average  costs  of  0.97  cents 
per  item  of  the  lot.  The  decision  losses  will  be  1.01  cents  per  item  of  the  lot  on 
the  average.  The  first  part  of  the  costs,  7.34,  depends  on  the  prior  distribution 
and  can  only  be  reduced  by  producing  (or  buying)  lots  of  better  quality.  The  second 
part,  1.98,  depends  on  the  sampling  plan  used.  Since  we  have  here  used  the  optimum 
plan  any  change  in  sample  size  or  acceptance  number  will  result  in  increased  costs. 
The  average  costs  of  accepting  all  lots  without  inspection  is  10.06  cents  per  item. 

The  two  functions  kQ(p)  =  Kq(p)/N  and  k(p)  ~  K(p)/N  have  been  shown  in  Fig.  1  for 
the  example  above . 


14.  General  remarks. 

There  exists  already  a  great  body  of  theories  and  tables  for  constructing  single 
sampling  attribute  plans  based  on  two  specified  quality  levels  (p^,P2)  and  some  further 
requirements.  To  see  how  the  present  paper  fits  into  this  the  most  important  systems 
have  been  listed  below  by  stating  the  "further  requirements"  for  each  system: 

(a) .  Specification  of  the  producer's  and  the  consumer's  risks,  see  for  instance  Peach 
and  Littauer  [7j  and  Grubbs  [_8  ]  . 

(b)  .  Specification  of  the  consumer's  risk  and  minimization  of  the  average  amount  of 
inspection  for  lots  of  process  average  quality  (p^)  in  the  case  of  rectifying 
inspection,  see  Dodge  and  Romig  [  9~j  . 

(c)  .  Specification  of  the  consumer's  risk  and  minimization  of  the  average  costs  for 
lots  of  process  average  quality  (p^),  i.e.  a  generalization  of  the  Dodge-Romig  LTPD 
system  requirering  specification  of  one  cost  parameter,  see  Hald  [lo]  . 
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(d)  .  Specification  of  two  cost  parameters,  p  and  p  ,  and  a  weight,  w„,  and 

r  s  * 

minimi  r*ntion  of  the  average  costs,  as  for  instance  in  the  present  paper. 

It  follows  from  the  results  of  the  present  paper  that  from  an  economic  point  of 
view  it  ir  not  advisable  to  fix  the  consumer's  or  the  producer's  risk.  On  the 
contrary  the  producer's  and  the  consumer's  risks  should  both  tend  to  »ero  with 
increasing  lot  size.  This  theorem  is  valid  not  only  for  the  double  binomial  prior 
distribution  but  for  any  prior  distribution,  and  it  is  valid  not  only  for  the 
Bayes  solution  but  also  for  the  minimax  solution  (p^<  p^<  )  t*1C  on^  difference 

being  the  speed  of  the  convergence.  For  a  discrete  prior  distribution  the  risks 
tend  to  zero  inversely  proportional  tc  N,  see  (63)  and  (64)  .  These  considerations 
lead  to  the  result  that  if  one  wants  a  system  with  a  fixed  risk  then  the  risk  should 
be  fixed  to  50  per  cent  at  a  point  between  p^  and  p^ .  We  may  therefore  increase  the 
list  of  systems  of  sampling  plans  above  by  the  following  item: 

(e)  .  Minimization  of  average  costs  for  lots  of  process  average  quality  (p^)  under 
the  restriction  that  P(pQ)  =  1/2.  Such  a  system,  named  the  IQL  system  (Indifference 
Quality  Level)  has  been  discussed  by  Hald,  see  [6]  and  [lo] ,  and  will  be  further 
discussed  in  a  forthcoming  paper.  This  system  requires  the  specification  of  p^  and 
a  cost  parameter.  In  view  of  the  asymptotic  relation  (66)  between  c  and  n  it  is 
clear  that  pQ  should  be  determined  from  (52)  . 

The  simplest  possible  system  based  on  the  specification  of  two  risks  and  having 
the  same  properties  as  the  Bayes  solution  may  be  formulated  as  follows: 

(f )  .  Specification  of  the  consumer's  or  the  producer's  risk  as  inversely  proportional 
to  lot  size,  and  P(p^)  “  1/2. 

This  system  requires  only  the  specification  of  one  parameter  (besides  the  two  quality 

levels)  and  it  is  extremely  simple  to  handle  both  mathematically  and  numerically. 

This  is  due  to  the  fact  that  the  equation  P(p  )  *  1/2  has  the  solution  c  =  np  +(p  -2)/3 

o  o  o 

(with  sufficient  accuracy  for  all  practical  applications,  perhaps  apart  from  the  case 
c  =  0  where  the  exact  solution  may  be  easily  found)and  that  the  other  equation, 

Q<P1)*S  a/N  say,  may  be  solved  with  respect  to  N  for  related  values  of  (n,c)  from 
the  first  equation.  Setting  c  =  0,1,2,...  and  solving  the  first  equation  for  n, 
the  second  equation  gives  N  =  a/(l-B(c,n, p^))  which  may  easily  be  found  by  means 
of  a  table  of  the  binomial  (or  the  Poisson)  distribution.  The  only  difficulty  lies 
in  the  choice  of  a.  If  the  problem  is  fully  specified  one  may  naturally  choose  a 
as  the  coefficient  of  1/N  in  (63)  and  the  system  will  then  asymptotically  give 
an  approximation  to  the  Bayes  solution.  The  reason  for  using  the  simple  system  will. 
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however,  usually  be  that  some  of  the  parameters  in.  the  problem  are  unknown  are  *  n 
that  case  the  choice  of  a  will  to  some  extort  be  arbitrary,  just  as  in  the  other  casc^ 
the  choice  of  the  producer's  or  the  consumer’s  risk  is  arbitrary.  This  system  ji 
sampling  plans  will  be  discussed  in  more  detail  in  the  forthcoming  paper  or.  th  IQI 
system. 

Turning  to  applications  it  is  important  to  notice  chat  a  system  of  sampling  plans 
in  practice  often  is  required  to  serve  several  purposes.  In  particular  we  shall 
here  stress  (a)  that  the  system  should  protect  the  consumer  against  deterioration 
of  the  prior  distributico,  (b)  that  the  system  should  work  as  an  incentive  for 
the  producer  to  produce  better  quality  or  at  least  to  keep  to  the  quality  agrecu 
upon,  see  Hill  [ll],  and  (c)  that  (average)  costs  should  be  minimized.  The  first 
two  requirements  are  concerned  with  consequences  of  changes  of  the  prior  distri¬ 
bution  and  the  problem  should  therefore  really  be  formulated  as  a  dynamic  cnc. 
However,  since  a  dynamic  model  at  present  is  lacking  wc  shall  try  to  indicate  hov 
the  Bayesian  solution  may  be  modified  to  take  requirements  (a)  and  (b)  into  account. 

One  of  the  arguments  advanced  against  the  Bayesian  method  in  general  has  been  that 
a  prior  distribution  does  normally  not  exist.  This  may  be  true  in  many  fields  but 
certainly  not  for  industrial  mass  production  with  its  effective  planning  and  control 
of  operations.  Admittedly  the  prior  distribution  may  change,  but  changes  are  usually 
rather  small  and  slow  within  a  given  production  period  in  which  the  same  machinery, 
techniques,  and  raw  materials  are  being  used.  We  arc  here  not  concerned  about  isolated 
very  poor  lots  which  may  occasionally  occur  since  any  sampling  plan  will  detect  such 
lots . 

Published  data  on  prior  distributions  arc  scarce.  Whether  the  double  binomial  distri¬ 
bution  is  a  reasonable  approximation  to  distributions  occurring  in  practice  is  not 
known.  According  to  the  experience  of  the  author  mixed  binomial  distributions  with 
beta-distributions  as  weight  functions  are  rather  common.  (A  paper  analogous  to  the 
present  one  will  present  the  corresponding  theory  and  tables  for  the  beta-distri¬ 
bution)  . 

One  of  the  drawbacks  of  the  Bayesian  solution  from  a  practical  point  of  viev;  is  that 
the  solution  may  be  acceptance  (or  rejection)  without  inspection.  If  one  is  rot 
completely  confident  that  the  prior  distribution  used  is  the  right  one  and  is 
then  a  sampling  plan  is  required  to  guard  against  deterioration  of  the  prio>*  distri¬ 
bution.  One  possibility  is  to  use  the  first  or  on<~  of  the  first  Bayesian  sampling 
plans  in  the  appropriate  table.  If  that  is  not  satisfactory  one  may  in  such  cases 
use  an  IQL  plan. 
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The  some  procedure  may  be  used  to  satisfy  requirement  (a)  above.  It  should  first 
of  all  be  noted  that  If  a  Bayesian  sampling  plan  exists  then  some  protection  against 
deterioration  of  the  prior  distribution  Is  automatically  obtained  and  the  protection 
may  In  the  usual  way  be  expressed  by  means  of  the  OC  curve.  It  is  always  easy  when 
the  plan  has  been  found  to  compute  the  consumer's  risk  and  then  to  decide  whether 
the  risk  Is  sufficiently  small.  If  the  consumer's  risk  is  too  large  one  may  again 
find  a  sufficiently  large  sample  in  the  same  table  or  turn  to  an  IQL  plan  or  a  LTPD 
plan. 

The  price  to  be  paid  for  obtaining  the  required  protection  is  naturally  that  the  plan 
used  will  not  minimize  costs  if  the  prior  distribution  holds.  If  the  change  in  the 
value  of  c  is  not  large  the  increase  in  costs  will,  however,  be  small. 

For  large  lots  the  consumer's  risk  for  the  Bayesian  sampling  plan  will  usually  be  much 
smaller  than  10  per  cent  so  that  the  problem  does  only  exist  for  small  lots. 

The  Incentive  for  the  producer  to  keep  to  the  specified  quality  is  usually  obtained 
by  alternating  between  normal  and  tightened  Inspection  in  a  specific  way  such  that 
the  system  reacts  upon  observed  changes  in  the  prior  distribution.  If  it  was  possible 
to  estimate  in  what  way  the  distribution  had  changed  the  reaction  could  be  made  to 
depend  on  the  change.  In  practice,  however,  one  wants  to  install  tightened  Inspection 
as  soon  as  possible  on  the  basis  of  some  over-all  criterion,  for  example  when  the 
number  of  lots  rejected  exceeds  some  critical  limit,  A  thorough  theory  does  not  exist 
but  some  rules  have  been  found  to  work  satisfactory  in  practice.  The  Military  Standard 
105  D  uses  the  same  sample  size  for  normal  and  tightened  inspection  and  a  reduced 
acceptance  number,  c^,  for  tightened  inspection.  The  difference  between  the  two 
acceptance  numbers,  c^-c^,,  equals  1  for  2  £  c^£  4,  2  for  5  £  c^S  20,  and  3  for  c^  21. 
For  c^=  0  or  1,  c^,  is  usually  equal  to  c^  but  the  sample  size  is  increased  for  tightened 
Inspection.  Similar  rules  may  be  used  for  the  present  tables  although  it  has  to  be 
realized  that  the  resulting  plans  will  not  be  minimum-cost  plans.  The  main  point  is, 
that  under  normal  conditions  the  plans  will  minimize  costs  and  that  the  plans  may  be 
adjusted  to  changes  in  the  prior  distribution  so  that  costs  are  minimized  under  the 
new  conditions.  If,  however,  the  incentive  aspect  of  sampling  inspection  is  more 
important  for  the  user  of  the  system  than  to  minimize  costs  in  case  of  change  of  the 
distribution  then  some  form  of  tightened  inspection  may  be  introduced  with  the  result 
that  during  periods  of  tightened  inspection  the  plans  will  not  minimize  costs. 
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Appendix, 


Master  Tables  of  Sampling  Plans. 

Tables  of  Conversion  Factors. 
Summary  of  Conversion  Formulas, 


All  sampling  plans  in  the  master  tables  assume  p^  ■  pg  and  =  0.05. 

In  the  first  set  of  tables  pf  *  10  X  and  (p^^)  take  on  the  values 

Pt  -  2.0,  2.5,  3.0,  3.5,  4.0,  5.0,  6.0,  7.0  7. 

P2  -  15.0,  17.5,  20.0,  25.0,  30.0  7,. 

In  the  second  set  of  tables  pf  ■  1  7.  and  (p^^)  take  on  tb«  values 

pt  -  0.20,  0.25,  0.30,  0.35,  0.40,  0.50,  0.60,  0.70  % 

p2  -  1.50,  1.75,  2.00,  2.50,  3.00,  3.50,  4.00,  5.00,  6.00,  7.00  7.. 


Single  Sampling  Tables  for 
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For  N  between  two  intervals  adjacent  in  the  table  find  (n,c)  for 
the  first  of  these  intervals  and  use  (n+1,c)  as  optimum  plan. 
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For  N  between  two  intervals  adjacent  in  the  table  find  (n,c)  for 
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Relation  between  and  *ox  fixed  (p^G#  p2c>  7£)» 


Use  the  sane  sampling  plan  for  V£  ■*  0.05  and  pfo  ■  0.01  (0.10)  as  for  V2  and 


(0. lOf)  where  f  is  given  in  the  table.  . 

?2  *  P2o/Pro  pl 
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Table  of  b^,  b^.*  and  b^, 


p2/pr 


1-4  n  n 


O'  Cn  ■£*  Ui  N>  sj  O'  m  MJ  N  M  O'  Ul  ^  U  M  >1  O'  Ui  ^  W  N 


P_  for  fixed  (p,,p,jV9) 


Use  N#  ■  Nf2  as  argument  in  the  master  tabic  to  find  (n*,c*). 

Pr  t:  pg  *  0.01X  or  0. 10X,  (p^P2^w2)  arc  Siven  in  the  master  tables, 
Pj^  =  lCXJpj  or  lOp^,  p2  =  100p2  or  10p2. 
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f 

') 


Correction  to  n*  due  to  a  change  in  p  =»  p_, 
z  •—  r  3 

Reference  value  p  «  p  -*  0.010,  n  =*  n*  +  g„. 

r  8  z 


P1 

X  -  0. 50 

0.60 

0.70 

0. 80 

0.90 

1.00 

1.25 

1.50 

1.75 

2.00 

3. 

0.2 

130 

100 

70 

50 

25 

0 

. 

0.3 

160 

120 

05 

55 

30 

0 

- 

- 

- 

- 

- 

0.4 

- 

150 

105 

65 

35 

0 

- 

- 

- 

- 

- 

0.5 

- 

- 

135 

85 

40 

0 

m 

- 

- 

- 
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0.6 

- 

- 

- 

115 

50 
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- 

- 

- 

- 
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- 

- 

- 

- 
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- 

- 

- 

- 

•  - 

0.2 

75 

55 

40 

25 

15 
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-30 

. 

_ 

- 

. 

0.3 

95 

70 

50 

30 

15 
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-35 

- 

- 

- 

- 

0.4 

- 

90 

60 

35 

15 
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-40 
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- 

- 
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0.5 

- 
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80 

45 

20 

0 

-45 
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- 
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25 
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-105 
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- 
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0.7 

m 
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40 
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-70 

-125 

- 

- 
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40 

30 

20 

15 

5 
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-25 

. 
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0.3 

55 

40 

25 

15 

5 

0 

-15 

-30 

-45 

- 

- 

0.4 

- 

50 

30 

20 

10 

0 

-20 

-35 

-50 

-65 

- 

0.5 

m 

- 

- 

25 

10 

0 

-20 

-40 

-55 

-70 

- 

0.6 

- 

- 

- 

- 

15 

0 

-25 

-45 

-60 

-80 

- 

0.7 

« 

- 

- 

- 

- 

0 

-30 

-55 

-75 

-90 

- 

0.2 

20 

15 

10 

5 

5 

0 

-  5 

-15 

-20 

-20 

_ 

0.3 

- 

20 

15 

10 

5 

0 

-10 

-15 

-20 

-25 

- 

0.4 

m 

- 

15 

10 

5 

0 

-1C 

-15 

20 

-25 

-45 

0.5 

- 

- 

10 

5 

0 

-10 

-20 

-25 

-30 

-50 

0.6 

- 

- 

- 

- 

- 

c 

-10 

-20 

-30 

-35 

-55 

0.7 

- 

- 

- 

- 

0 

-15 

-25 

-35 

-40 

-60 

0.2 

15 

10 

5 

5 

0 

0 

-  5 

-10 

-10 

-15 

-25 

0.3 

- 

15 

10 

5 

0 

0 

-  5 

-10 

-10 

-15 

-25 

0.4 

- 

- 

- 

5 

5 

0 

-  5 

-10 

-15 

-15 

-25 

0.5 

- 

- 

- 

5 

0 

-  5 

-10 

-15 

-20 

-30 

0.6 

•• 

- 

- 

- 

- 

0 

-10 

-15 

-20 

-20 

-35 

0.7 

' 

0 

-10 

-15 

-20 

-25 

-35 

For  p  * 
rr 

•o 

CO 

II 

0.10  the  correction  is  ^ 

/10 

(rounded  down) 

• 

1.5 


2.0 


3.0 


5.0 


7.0 


Conversion  factor  for  N  due  to  a  change  in  w2 


Reference  value  of  w0  =  0.05,  p  -  p  . 

2  s  r 

Use  N*  *  Nf^  as  argument  in  the  master  tabic  to  find  (n*,c*)„ 

P2/Pr 


100w2 

1.5 

2.0 

3.0 

5.0 

7.0 

1 

J 

1 

0. 54 

0.  56 

0.  58 

0.61 

0.63 

O* 

« 

0 

0.46 

0.48 

0.51 

0. 54 

0.57 

0  7 

2  i 

0.70 

0.72 

0.  74 

0.76 

0.77 

0,2 

i 

0.65 

0, 66 

0.68 

0.71 

0.73 

0.7 

j 

**  ! 

0.82 

0.83 

0.84 

0. 06 

0.07 

0.2 

j 

0.78 

0.79 

0.81 

0.83 

0. 84 

0.7 

L  1 

0.92 

0. 92 

0.  93 

0.94 

0. 94 

0.2 

4 

} 

0. 90 

0. 90 

0.91 

0.92 

0.  93 

0.7 

1.00 

1.00 

1.00 

1.00 

0 

0 

• 

0.2 

j 

1.00 

1.00 

1.00 

1.00 

1.00 

0  7 

1.07 

1.07 

1.06 

1.06 

1.05 

0.2 

0  s 

1.09 

1.09 

1.08 

1.07 

1.06 

0.7 

7 

1.14 

1.13 

1.12 

1.10 

1.09 

0.2 

1.17 

1.16 

1.15 

1.13 

1.12 

0.7 

8 

1.19 

1.18 

1.16 

1.15 

1.13 

0.2 

1.25 

1.24 

1.21 

1.19 

1.17 

0.7 

Q 

1.25 

1.23 

1.21 

1.18 

1.17 

0.2 

y 

1.32 

1.30 

1.27 

1.24 

1.22 

0.7 

10 

1.30 

1.27 

1.25 

1.22 

1.2C 

0.2 

1.39 

1.37 

1.33 

1.29 

1.26 

0.7 

12 

1.38 

1.36 

1.32 

1.28 

1.25 

0.2 

1.51 

1.48 

1.43 

1.38 

1.34 

0,7 

14 

1.46 

1.43 

1.38 

1.33 

1.30 

0.2 

1.63 

1.58 

1.52 

1.45 

1.40 

0.7 

16  ! 

1.53 

1.49 

1.44 

1.38 

1.34 

0.2 

1.73 

1.68 

1.61 

1.52 

1.46 

0.7 

18  | 

1 

1.60 

1.55 

1.49 

1.42 

1.38 

0.2 

1.83 

1.77 

1.68 

1.58 

1.52 

0.7 

l 

20  | 

i 

1.65 

1.60 

1.53 

1.46 

1.41 

0.2 

1.92 

1.85 

1.75 

1.64 

1.  56 

0,7 
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Correction  ^  to  n*  due  to  a  change  in  ”2* 

Reference  value  of  v,  -  0.05,  p  ■  p  ■  0,01,  n  ■  n*  +  g, . 

z  9  r  * 

p2/pr 


100w2 

1.5 

2.0 

3.0 

5.0 

7.0 

pl/pr 

1 

-123 

-90 

-60 

-35 

-25 

0.2 

1 

-205 

-125 

-70 

-35 

-25 

0.7 

O 

-70 

-50 

-35 

-20 

-15 

0.2 

4 

-115 

-70 

-40 

-20 

-15 

0.7 

-40 

-30 

-20 

-10 

-10 

0.2 

J 

-65 

-40 

-25 

-10 

-10 

0.7 

A 

-20 

-15 

-10 

-5 

-5 

0.2 

-30 

-20 

-10 

-5 

-5 

0.7 

c 

0 

0 

0 

0 

0 

0.2 

0 

0 

0 

0 

0 

0.7 

4 

15 

10 

5 

5 

5 

0.2 

25 

15 

10 

5 

5 

0.7 

7 

25 

20 

15 

5 

5 

0.2 

i 

45 

25 

15 

10 

5 

0.7 

8 

40 

30 

20 

10 

5 

0.2 

60 

40 

20 

10 

10 

0.7 

q 

50 

35 

20 

15 

10 

0.2 

-r 

80 

50 

25 

15 

10 

0.7 

10 

55 

40 

25 

15 

10 

0.2 

90 

55 

30 

15 

10 

0.7 

12 

75 

50 

35 

20 

15 

0.2 

120 

70 

40 

20 

15 

0.7 

14 

85 

60 

40 

25 

15 

0.2 

140 

85 

50 

25 

15 

0.7 

16 

100 

70 

45 

25 

20 

0.2 

160 

100 

55 

30 

20 

0.7 

18 

110 

80 

50 

30 

20 

0.2 

175 

110 

60 

30 

20 

0.7 

20 

120 

85 

55 

30 

20 

0.2 

195 

120 

65 

35 

25 

0.7 

For  p  ■  p  *  0.10  the  correction  is  g./10  (rounded  down). 
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Summary  of  conversion  formulas 
to  find 


(n,c)  corresponding  to  (N,  pf,  pg,  p^  P2, 
from 

(n*,c#)  in  the  master  table  for  (N*,  pro, 


For 


> 


0.05  'i 
0.05  j 


use  master  table  with  p 

rro 


0.01 

0.10 


X 

s 


(i— &.Y1 

V  VV*5!*  J 


Formula  1. 

w7(P2-Pr)  P2+19r2P, 

7,  «  — ; - r  and  Xp  *  ■  •■rrg -  • 

2  wi(Pr‘Pi>  ro  1+1 972 

Use 

n*  -  m.a\,  pro,  p1o  -  Pj/x,  p2o-  p2/x 

as  arguments  to  find  (n*,c*)  in  the  master  table. 

(n,c)  -  (n*/X,c*). 

If  (Pj0^P2O)  outside  the  tabulated  range  use  formula  2. 

Formula  2. 

X  -  pr/pr0,  Cj  -  Pj/pr,  P2  -  P2/Pr. 

Use 

N*  -  NX8Xf1(w2,p1,p2),  pro,  Pu  ■  Olpro,  p2o 


as  arguments  to  find  (n*,c*)  in  the  master  table 
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